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Integrated density of states and Wegner estimates for 
random Schrodinger Operators 

Ivan Veselic 



Abstract. We survey recent results on spectral properties of random Schro- 
dinger operators. The focus is set on the integrated density of states (IDS). 
First we present a proof of the existence of a self-averaging IDS which is general 
enough to be applicable to random Schrodinger and Laplace-Beltrami opera- 
tors on manifolds. Subsequently we study more specific models in Euclidean 
space, namely of alloy type, and concentrate on the regularity properties of 
the IDS. We discuss the role of the integrated density of states and its reg- 
ularity properties for the spectral analysis of random Schrodinger operators, 
particularly in relation to localisation. Proofs of the central results are given 
in detail. Whenever there are alternative proofs, the different approaches are 
compared. 

Resumen. Revisamos resultados recientes en propiedades espectrales de oper- 
adores de Schrodinger aleatorios. Nos enfocamos principialmente en la densi- 
dad integrada de estados (IDS). Primero presentamos una prueba de la existen- 
cia de la IDS y su propriedad auto-promediadora (self-averaging). El metodo 
es suficientemente general para ser aplicable a operadores de Schrodinger y de 
Laplace-Beltrami aleatorios en varicdades de Riemann. Posteriormente estu- 
diamos los modelos mas especificos en el espacio Euclidiano, a saber de tipo 
aleacion, y nos concentramos en las propiedades de la regularidad de la IDS. 
Discutimos el papel de la densidad integrada de estados y sus propiedades de 
regularidad para el analisis espectral de operadores de Schrodinger aleatorios, 
particularmente con relation a localization espectral. Las pruebas de los re- 
sultados centrales son descritas en detalle. Cuando hay pruebas alternativas, 
los enfoques diferentes se comparan. 
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1. Random operators 

1.1. Physical background. Random Schrodinger operators are used as mod- 
els of disordered solids within the framework of quantum mechanics. 

A macroscopic solid consists of an order of magnitude of 10 23 of nuclei and 
electrons. The resulting Hamiltonian taking into account all interactions is highly 
complicated. To arrive at a Schrodinger operator which can be studied in some 
detail one neglects the electron-electron interaction and treats the nuclei in the 
infinite mass approximation. Thus one arrives at an one-electron Schrodinger oper- 
ator with an external potential due to the electric forces between the electron and 
the nuclei, which are assumed to be fixed in space. 

In the case that the nuclei are arranged periodically on a lattice, the potential 
energy of the electron is a periodic function of the space variable. 
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On the other hand, the electron could be moving in an amorphous medium, 
in which case there is no large group of symmetries of the Hamiltonian. However, 
from the physical point of view it is reasonable to assume that the local structure of 
the medium will be translation invariant on average. This means that we consider 
the potential which the electron experiences as a particular realisation of a random 
process and assume stationarity with respect to some group of translations. More- 
over, physical intuition suggests to assume that the local properties of the medium 
in two regions far apart (on the microscopic scale) are approximately independent 
from each other. Therefore the stochastic process describing the potential should 
have a correlation function which decays to zero, or — more precisely — should be 
crgodic. 

There are interesting models which lie between the two extreme cases of lattice- 
periodic and amorphous media. They still have an underlying lattice structure 
which is, however, modified by disorder. Probably the best studied Hamiltonian 
with this properties is the alloy type model. We leave its precise definition for the 
next paragraph and introduce here a special case on the intuitive level. Consider 
first the potential 

Vu(x) := ^ u k(v,%) 

k& d 

Each k corresponds to a nucleus sitting on a lattice point. The function Uk(u>, •) 
describes the atomic or nuclear potential at the site k and depends on the random 
parameter lo which models the different realisations of the configuration of the 
nuclei. If there is only one type of atom present, which has a spherically symmetric 
potential, all the •) are the same, and V u is periodic. Now assume that there 

are two kinds a and b of atoms present, which have spherically symmetric atomic 
potentials of the same shape, but which differ in their nuclear charge numbers. 
In this case the potential looks like 

V u (x) := q a u(x~k) + q b u(x - k) 

k occupied by a k occupied by b 

If the two sorts of atoms are arranged on the lattice in a regular pattern, this again 
gives rise to a periodic potential. 

However, there are physically interesting situations (e.g. binary alloys) where 
the type of atom sitting on site k is random, for example obeying the law 

P{fc is occupied by atom a} — P, P{k is occupied by atom 6} = 1 — P P > 

Here P{. . . } denotes the probability of an event. If we furthermore assume that the 
above probabilities are independent at each site and the parameter P is the same 
for all k, we arrive at the continuum Bernoulli-Anderson potential 

V ul (x) = ^ 9k M u{x - k) 

k 

Here qk(u>) € {q a ,qb},k € Z d denotes a collection of independent, identically dis- 
tributed Bernoulli random variables and u is some atomic potential. 

This model is a prototype which has motivated much research in the physics 
and mathematics literature, a part of which we will review in the present work. 

Properties of disordered systems are discussed in the books |35L 1911 1208) from 
the point of view of theoretical physics. The mathematical literature on random 
Schrodinger operators includes the books [471 1241) , the introductory article (158] , 
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a section on random Jacobi matrices in |65| . the Lifshitz memorial issue |207| . and 
a monograph specialised on localisation phenomena 295 . 

In the present work we focus on the spectral features of random Hamiltonians 
which are encoded in, or in close relation to the properties of the integrated density 
of states (IDS). We present the proofs of most of the results in detail. Indeed, for 
the two central theorems — the existence of an self-averaging integrated density 
of states and for the Wegner estimate — two independent proofs are given and 
compared. 



1.2. Model and Notation. Let us start with some mathematical notation. 
The symbols R, Z, N, No denote the set of reals, the set of integers, the set of natural 
numbers, and the set on non-negative integers, respectively. For a set A we denote 
by A c its complement. An open subset of R d will be denoted by A, and if there is a 

sequence of such sets its members will be denoted Ai, . . . , A;, The symbol |A| 

is used for the Lebesgue measure of A. We write |x| for the norm ofxeR d , while 
the norm of a vector / in a function space is denoted by ||/||. 

The Hilbert space of (equivalence classes of) measurable functions on A which 
are square integrable with respect to Lebesgue measure is denoted by L 2 (A). Simi- 
larly, L P (A) with p > stands for the space of measurable functions / such that \f\ p 
is integrable, while L°°(A) is the set of measurable functions which are essentially 
bounded with respect to Lebesgue measure. The space of sequences {a„}„ e N such 
that |a n | p is summable is denoted by l p (N). Note that the case p £]0, 1[ is included 
in our notation. In our context we will often choose the exponent p dependent on 
the dimension of the configuration space. In the following we denote by p(d) any 
number in [1, oo[ which satisfies 



(1.1) p{d) 



> 2 if d < 3, 

> d/2 if d > 4 



The symbols C(A), C°°(A) stand for the continuous, respectively smooth, functions 
on A. The subscript c in C C (A), C£°(A), L P (A) means that we consider only those 
functions which have compact support in A. In the sequel we will often consider 
potentials from the class of functions which are uniformly locally in L p . More 
precisely, / is in the set of uniformly locally L p -functions, denoted by £u nif i oc (R d ), 
if and only if there is a constant C such that for each y g R d 

\f(x)\ p dx<C 



1, 



\x-y\<l 

The infimum over all such constants C is ||/||p unif loc . 

Let A denote the Laplacian on R d . If we choose its operator domain T>(A) to 
be the Sobolev space Ly|(R d ) of functions in L 2 (R d ) whose second derivatives (in 
the sense of distributions) are square integrable, it becomes a selfadjoint operator. 
The restriction of A to a true open subset A C R d becomes selfadjoint only if we 
specify appropriate boundary conditions (b.c). Dirichlet b.c. are defined in Remark 
12.2.31 For the definition of Neumann and periodic b.c. see for instance |256| . 

Let A, B be two symmetric operators on a Hilbert space 7i, whose norm we 
denote by || ■ ||. We say that B is (relatively) A-bounded if the domains obey the 
inclusion V(A) C T>(B) and there are finite constants a and c a such that for all 
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/ G V{A) 

(1-2) ||5/||<a||A/|| + c ||/|| 

The infimum over all a such that the estimate holds with some c a is called relative 
bound (of B with respect to A). If B is A-bounded with relative bound zero, 
we call it infinitesimally A-bounded. Let A be selfadjoint, and B symmetric and 
relatively A-bounded with relative bound smaller than one. Then the operator 
sum A + B on the domain T>(A) is selfadjoint by the Kato-Rellich Theorem. We 
will apply this result to the sum of the negative Laplacian and a potential. A 
multiplication operator by a function V £ ifLjf i oc (R d ) is infinitesimally A-boundcd 
if p = p(d), cf. [2561 Thm. XIII. 96]. Moreover, the constant c a in fl~2|l depends 
only on ||y||p,unif,ioc- Thus the sum H := -A + V is selfadjoint on W 2 2 (R d ). If B 
is relatively A-boundcd in operator sense with relative bound a it implies 

(/, Bf) < a(f, Af) + C a (f, /) for some C a £ R 

which is called relative form-boundedness of the form Qb (/,<?) := (f,Bg) with 
respect to Qa{L g) '■= (f,Ag). See § VI. 1.7 in |152| for more details. 

The triple (O, Bn,¥) stands for a probability space with associated cr-algebra 
and probability measure, while E{. ..} denotes the expectation value with respect 
to P. A collection Tj : f2 — > fi, j £ J of measure preserving transformations is 
called ergodic if all measurable sets in fl which are invariant under the action of all 
Tj,j £ J have measure zero or one. 

Definition 1.2.1. Let p = p(d) be as in i|l.l|) . u £ LP(R d ) and qt.: — > 
R, k £ Tj d be a sequence of bounded, independent, identically distributed random 
variables, called coupling constants. Then the family of multiplication operators 
given by the stochastic process 

(1.3) V u {x) := ^2 qk{u)u(x-k) 

kez d 

is called alloy type potential. The function u is called single site potential. Let 
H := — A + V pcl be a periodic Schrodinger operator with V peT £ L^ ni{ i oc (R d ). The 
family of operators 

(1.4) H U :=H + V U , oj£il 

is called alloy type model. 

The distribution measure of the random variable go will be called single site 
distribution and denoted by /i. If not stated otherwise, in the sequel we assume 
that p is absolutely continuous with respect to the Lebesgue measure and has a 
bounded density. The density function is denoted by /. 

Due to our assumptions on the boundedness of the coupling constants, for each 
a > there is a constant c a such that for all w and all ip £ T>(A) 

HKVII < all A^|| + C»||^||, HVperV'll < o||A^|| + C»|W 

In particular Hq and all are selfadjoint on the operator domain of A. It will 
be of importance to us that the constant c a may be chosen independently of the 
random parameter u>. 
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Remark 1.2.2. (a) In several paragraphs we study Hamiltonians as in Defi- 
nition 11.2.11 but where some of the hypotheses on the single site potential or the 
coupling constants are relaxed. More precisely, we will consider single site po- 
tentials with non-compact support and coupling constants which are unbounded, 
correlated, or do not have a bounded density. 

(b) If the coupling constants are not bounded, one has to impose some moment 
condition to make sure that the alloy type model still makes sense. The main 
difference (to the bounded case) is that for w in a set fl' C SI of full measure the 
operator H u will be (essentially) selfadjoint, however this will fail to hold for w in 
the complement SI \ Of . See for example |162L 11631 1165| for more details. 

(c) There is a group of measure preserving transformations T^,k £ Z d on 
(f2,Bn,P) such that Jd obeys 

V u (x-k)=V 7ilU {x) 

In other words the stochastic process V : SI x R — > R is stationary with respect to 
translations by vectors in 1 d . Moreover, the group T k , k £ 1* d acts ergodically on 
SI, therefore we call V an *Z d -ergodic potential. 

To see that the above statements are true we pass over to the canonical proba- 
bility space SI = x fc6Z d R, equipped with the product measure P := ®kei d t 1 - Now 
the stochastic process {irk}kez d i defined by nk(uj) = uik for all k £ Z d , has the same 
finite dimensional distributions as {qk}k- It is easily seen that the transformations 
(Tk(u)))j :— ujj-k are measure preserving and that the group Z d acts ergodically on 

si. 

Using the stochastic process {~Kk\k the alloy type potential can be written as 

(1.5) Vu(x) := u k u(x- k) 

kez d 

which we will use without distinction in the sequel. 

Abstracting the properties of stationarity and ergodicity we define general ran- 
dom potentials and operators with Z d -ergodic structure. 

Definition 1.2.3. Let V: SI X R d — * R be a stochastic process such that for 
almost all uj £ SI the realisation of the potential obeys V u £ L^ nit loc (R d ), p = p(d) 
and additionally E{||VL,XA||p} < oo, where A is an unit cube. Let T^,k £ II 1 be a 
group of measure preserving transformations acting ergodically on (Sl,Z?n,P) such 
that 

V u {x-k) =V TkU (x) 

Then we call {V^} w a {X d -ergodic) random potential and {H^}^ with H w = — A+VL 
a (1i d -ergodic) random operator. 

The restriction of H w to an open subset A will be denoted by if we impose 
Dirichlet boundary conditions and by H^' N in the case of Neumann b.c. While we 
will be mainly concerned with Z d -ergodic operators we will give some comments 
as asides on their counterparts which are ergodic with respect to the group R d . 
The recent overview 203 is devoted to such models that model amorphous media. 
Insight in the research on almost periodic operators can be obtained for instance 
in |298l 12991 EB E3 1241] and the references therein. 

Remark 1.2.4. All Z d -ergodic potentials can be represented in a form which 
resembles alloy type potentials. Namely, for such V: SI x R d — > R there exists a 
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sequence fu , k € % d of random variables on O taking values in the separable Banach 
space L p (R d ) such that V can be written as 



This representation is of interest because it ensures that after passing to an equiv- 
alent probability space and stochastic process one may assume that the sigma al- 
gebra on £1 is countably generated. See [156] and Remark 2.8 in j!99| for more 
information. 

1.3. Transport properties and spectral types. The main interest in the 
study of random operators is to understand the transport properties of the media 
they model. In the particular case of the quantum mechanical Hamiltonian of an 
electron in a disordered solid the electric conductance properties are the main object 
of interest. 

The Hamiltonian governs the equation of motion, i.e. the time dependent Schro- 
dinger equation 



The time evolution of the vector ip(t) in Hilbert space describes the movement of 
the electron. Since we chose the space representation in the Schrodinger picture, 
we can can think of ip(t) as a wave packet which evolves in time. The square 
of its absolute value \ip(t,-)\ 2 S L 1 (M d ) is a probability density. More precisely, 
J A \ip(t, x)\ 2 dx is the probability to find the electron in the set A C R d at time t. 

For a given initial state ■f/'o := ip{0) supported in a compact set A one would 
like to know whether for large times the function ip{t) stays (essentially) supported 
near A, or moves away to infinity. In the first case one speaks of a bound state, 
since it remains localised near its original support for all times. The other extreme 
case would be that ip(t) leaves any compact region in M. d (and never comes back) 
as time goes to infinity. Such a state is called a scattering or extended state. By 
the RAGE theorem, cf. e.g. [2571 1651 1295) . it is possible to relate the dynamical 
properties of states just described to the spectral properties of the Hamiltonian. 
Roughly speaking, bound states correspond to pure point spectrum and scattering 
states to (absolutely) continuous spectrum. For a more precise statement consult 
for instance (25T1 16311295] , 

This motivates the systematic study of spectral properties of the introduced 
Schrodinger operators. If a random Schrodinger operator exhibits almost surely 
only pure point spectrum in an energy region one speaks of Anderson or spectral 
localisation. The name goes back to Anderson's seminal work |13| . This property 
has been established for a variety of random models. In most of those cases one 
can additionally prove that the corresponding eigenfunctions decay exponentially in 
configuration space, a phenomenon called exponential (spectral) localisation. The 
situation is different for random potentials with long range correlations, where 
sometimes only power-law decay of the eigenfunctions has been established [1681 
H09ll329] . 

If an energy interval contains almost surely only pure point spectrum, we call it 
localisation interval. An eigenfunction of H u which decays exponentially is called an 
exponentially localised eigenstate. The region or point in space where the localised 



(1.6) 




104 



I. VESELIC 



state has its highest amplitude will be called localisation centre (we will not need 
a mathematically precise definition of this notion) . 

However, it turns out that the spectrum captures only a rough view on the 
dynamical properties of the quantum mechanical system. A more detailed un- 
derstanding can be obtained by studying the time evolution of the moments of 
the position operators. This led to a formulation of several criteria of dynami- 
cal localisation. The strongest characterisation of this phenomenon, namely strong 
dynamical localisation in Hilbert- Schmidt topology means that for all q > 



Here P^(I) denotes the spectral projection onto the energy interval / associated to 
the operator H\ }1 || • \ hs denotes the Hilbert-Schmidt norm, K C M. d is any compact 
set, and \X\ denotes the operator of multiplication with the function \x\. For the 
interpretation of l|1.7|) as non-spreading of wave-packets one chooses f(y) = e~ tty . 
Dynamical localisation (|1.7|) implies in particular that the random Hamiltonian 
H w exhibits spectral localisation in /. The reader can consult e.g. |262t 11181 12951 
11201 11191 Rj] to get an insight how the notion of dynamical localisation evolved and 
for recent developments. In |262l I70L I148j examples are discussed where spectral 
localisation occurs, but certain dynamical criteria for localisation are not satisfied. 

For the purpose of the present paper these distinctions are not crucial. In 
the case of the alloy type model, to which we devote most attention, spectral and 
dynamical localisation coincide, cf. [6711119) . In the sequel we mean by localisation 
that the considered operator exhibits in a certain energy interval only pure point 
spectrum, and that the corresponding eigenfunctions decay sufficiently fast. 

Since we are dealing not just with a single Hamiltonian, but with a whole family 
of them, we have to say something on how the spectral properties depend on the 
parameter u> describing the randomness: most properties of the spectrum of an 
operator pertaining to the family {H^,}^ hold almost surely, i.e. for w in a set such 
that its complement has measure zero in f2. This is at least true for the properties 
we discuss in the present paper. 

We shortly describe what kinds of spectral types one expects from the physical 
point of view for random Schrodinger operators, say of alloy type. In case there are 
rigorous results which have confirmed this intuition we quote the reference. 

In one space dimension the spectrum is pure point for all energies almost surely. 
Rigorous proofs of this statement can for instance be found in I128L I193j . 

In three or more dimensions it is expected that the spectrum is pure point 
near the boundaries of the spectrum while in the interior it is purely absolutely 
continuous. In the latter case one speaks also of an energy region with delocalised 
states. However, for alloy type Hamiltonians the proof of delocalisation is open. 
Some results on existence of absolutely continuous spectrum for random models of 
a different type can be found in [1761 11021 1160) . The two regions with localised, 
respectively delocalised states are separated by a threshold, the so called mobility 
edge, for partial results see [1471 11601 11191 1103} . 

The literature on the existence of pure point spectrum is extensive. We discuss 
it in more detail in S 13.21 

How large the intervals with point or continuous spectrum are, depends on the 
disorder present in the model. For instance, in (|1 -4|) one could introduce a global 
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coupling constant A in front of the potential 

Huj = Hq + XVu! 

Now large A means large disorder, small A small disorder. The larger the disorder, 
the larger is the portion of the spectrum which contains localised states. For other 
types of random Schrodinger operators there are similar ways to introduce a disorder 
parameter. 

The phenomenon that localised states emerge at the edges of the spectrum 
can be understood in terms of the so-called fluctuation boundaries. These are the 
regions of the spectrum which correspond to extremely rare configurations of the 
potential. Consequently, the density of states (or the spectral density function, see 
the next section for a precise definition) is very thin in this region. This has been 
first understood on physical grounds by Lifshitz. Today the tails of the density 
of states at the fluctuation boundaries bear the name of Lifshitz- asymptotics or 
Lifshitz-tails. We give a (non-exhaustive) list of works devoted the study of this 
asymptotics: [238l 12391 12351 ITl5l Il64l 12801 UtTf! 12221 12231 l28ll l3(75l 12241 
ITIT1 12931 H84l IT851 Il83l Il82] and the references at the end of § E3 

The existence of localised states for random Schrodinger operators is in sharp 
contrast to the features of periodic operators. Indeed, for operators with periodic 
potential, satisfying some mild regularity assumptions, it is known that the spec- 
trum is purely absolutely continuous, [321 12881 1273L 12 741 1195| . This difference 
might seem somewhat surprising, given the similarity of the structure of an alloy 
type and a periodic operator. 

1.4. Outline of the paper. In the present section we discussed the physical 
motivation to study random operators and introduced the models we will analyse 
in the sequel. The next section is devoted to the proof of the existence of the 
IDS for these models. In the third section we discuss why regularity properties of 
the IDS are of interest. The following two sections give two independent proofs 
of the continuity of the IDS. Both approaches are suitable to extensions in various 
directions. We review some recently obtained results. For more details see the table 
of contents. 
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the IIMAS. This is an opportunity to thank D. Hundcrtmak, R. Killip, V. Kostrykin, 
D. Lenz, N. Peyerimhoff, and O. Post for the pleasure of collaborating with them on 
some of the results presented here, and in particular W. Kirsch who introduced me 
to the topic of random Schrodinger operators. I am grateful to E. Giere, P. Miiller, 
R. Muno, K. Schnee, K. Veselic, S. Warzel and the referees for valuable comments. 
Financial support by a fellowship of the DFG and hospitality at CalTech by B. Si- 
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2. Existence of the integrated density of states 

Intuitively, the integrated density of states (IDS) measures how many electron 
energy levels can be found below a given energy per unit volume of a solid. It can 
be used to calculate the free energy and hence all basic thermodynamic quantities 
of the corresponding non-interacting many-particle system. 
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To define the IDS mathematically one uses an exhaustion procedure. More 
precisely, one takes an increasing sequence A; of open subsets of R d such that each A/ 
has finite volume and [j. A; = R d . Then the operator H^, which is the restriction of 
H^, to A; with Dirichlet boundary conditions, is selfadjoint, bounded below and its 
spectrum consists of discrete eigenvalues Ai(fl^) < ^{H 1 ^) < • • ■ < A„(i/^) — > oo. 
Here A„ = A„+i means that the eigenvalue is degenerate and we take this into 
account in the enumeration. 

The normalised eigenvalue counting function or finite volume integrated density 
of states N}, is defined as 

(2.1) K(E) := *{n\K(Hl)<E} 
\ ) u \ ) | A; | 

The numerator can equally well be expressed using the trace of a spectral projection 
#H K(H l J <E} = Tr[p< (] - oo,E[) 

Note that N l u : R ->■ [0, oo[ is a distribution function of a point measure for all I G N, 
i.e. N^E) — — oo,E[). Here ir u is the finite volume density of states measure 
defined by 

By definition a distribution function is non-negative, left-continuous and monotone 
increasing. In particular, it has at most countably many points of discontinuity. 

Under some additional conditions on the random operator and the exhaustion 
sequence A; , I £ N one can prove that 

(i) For almost all oj € f2 the sequence A^, converges to a distribution function 
N u as / goes to infinity. This means that we have N^(E) — > N U (E) for all 
continuity points E of the limit distribution N u . 

(ii) For almost all u> € fi the distribution functions N u coincide, i.e. there is 
an ^-independent distribution function N such that A^ = for almost all 
uj. This function A^ is called the integrated density of states. Note that its 
independence of uj is not due to an explicit integration over the probability 
space J7, but only to the exhaustion procedure. This is the reason why the 
IDS is called self-averaging. 

(iii) In most cases there is a formula for the IDS as an expectation value of a trace 
of a localised projection. For Z d -ergodic operators it reads 

(2.2) N(E) := E{Tr[xAP w (]-oo,£[)]} 

Here A denotes the unit box ]0, l[ d , which is the periodicity cell of the lattice 
Z d . Actually, one could choose certain other functions instead of xa, yielding 
all the same result, cf. Formula l|2.15(l . The equality (|2.2(l holds for R d -ergodic 
operators, too. It is sometimes called Pastur-Subin trace formula. 

In the following we prove the properties of the IDS just mentioned by two meth- 
ods. In 12.21 - -- 12.61 a complete proof is given using the Laplace transforms of the 
distribution functions N^, while § 12. 71 is devoted to a short sketch of an alternative 
method. It uses Dirichlet-Neumann bracketing estimates for Schrodinger opera- 
tors, which carry over to the corresponding eigenvalue counting functions. These 
are thus super- or subadditive processes to which an ergodic theorem [12| can be 
applied. 
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Actually the proof using Laplace transforms will apply to more general situ- 
ations than discussed so far, namely to more general geometries than Euclidean 
space. To be precise, we will consider random Schrodinger operators on Riemann- 
ian covering manifolds, where both the potential and the metric may depend on 
the randomness. This includes random Laplace-Beltrami operators. 

We follow the presentation and proofs in |243L [199 . The general strategy 
we use was developed by Pastur and Subin in |237| and |298j for random and 
almost-periodic operators in Euclidean space. A particular idea of this approach 
is to prove the convergence of the Laplace transforms C 1 ^ of the normalised finite 
volume eigenvalue counting functions instead of proving the convergence of Njj 
directly. This is actually the main difference to the second approach we outline in 
§ 12.71 which is taken from |161j . The Pastur-Subin strategy seems to be better 
suited for geometries with underlying group structure which is non-abelian. 

Indeed, one of the differences between random operators on manifolds and those 
on K d is that the operator is equivariant with respect to a group which does not 
need to be commutative. This means that one has to use a non-abelian ergodic 
theorem to derive the convergence of the distribution functions N l u or, alternatively, 
of their Laplace transforms C l u . 

This imposes some restriction on the strategy of the proof since the ergodic 
theorems which apply to non-abelian groups need more restrictive assumptions 
than their counterparts for commutative groups, cf. also Remark 12. 6. 21 

For processes which are not additive, but only super- or subadditive, there is a 
non-abelian maximal ergodic theorem at disposal (cf. 6.4.1 Theorem in |194j ) but 
so far no pointwise theorem. This is also the reason why the Dirichlet-Neumann 
bracketing approach of § 12.71 does not seem applicable to random operators liv- 
ing on a covering manifold with non-abelian deck-transformation group (covering 
transformation group). 

2.1. Schrodinger operators on manifolds: motivation. In this section 
we study the IDS of random Schrodinger operators on manifolds. Let us first 
explain the physical motivation for this task. 

Consider a particle or a system of particles which are constrained to a sub- 
manifold of the ambient (configuration) space. The classical and quantum Hamil- 
tonians for such systems have been studied e.g. in |226UllU] (see also the references 
therein). To arrive at an effective Hamiltonian describing the constrained motion 
on the sub-manifold, a limiting procedure is used: a (sequence of) confining high- 
barrier potential(s) is added to the Hamiltonian defined on the ambient space to 
restrict the particle (system) to the sub-manifold. In [2261 1110] one can find a dis- 
cussion of the similarities and differences between the obtained effective quantum 
Hamiltonian and its classical analogue. 

A important feature of the effective quantum Hamiltonian is the appearance of 
a so-called extra-potential depending on the extrinsic curvature of the sub-manifold 
and the curvature of the ambient space. This means that even if we disregard 
external electric forces the relevant quantum mechanical Hamiltonian of the con- 
strained system is not the pure Laplacian but contains (in general) a potential 
energy term. This fact explains the existence of curvature-induced bound states in 
quantum waveguides and layers, see jlOOL 1891 12111 190) and the references therein. 

As is mentioned in 226 , the study of effective Hamiltonians of constrained 
systems is motivated by specific physical applications. They include stiff molecular 
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bonds in (clusters of) rigid molecules and molecular systems evolving along reaction 
paths. From the point of view of the present paper quantum wires, wave guides and 
layers are particularly interesting physical examples. Indeed, for these models (in 
contrast to quantum dots) at least one dimension of the constraint sub-manifold is of 
macroscopic size. Moreover, it is natural to assume that the resulting Hamiltonian 
exhibits some form of translation invariance in the macroscopic direction. E.g. it 
may be periodic, quasiperiodic or — in the case of a random model — ergodic. 

For random quantum waveguides and layers the existence of dense point spec- 
trum is expected, cf. the discussion of localisation in Paragraph F° r a waveg- 
uide embedded in Euclidean space this has been rigorously proven in |174l I175| . 
The question of spectral localisation due to random geometries has been raised al- 
ready in ffifl] . There the behaviour of Laplace-Beltrami operators under non-smooth 
perturbations of the metric is studied. 

A second motivation for the analysis of operators on manifolds studied in 
this section comes from differential geometry. The spectral properties of periodic 
Schrodinger operators on manifolds have attracted the interest of various authors. 
A non-exhaustive list is [M l30ll Il86l 13021 13031 HT1 Hoi 13001 Il50l 124411245] , 
By 'periodic' we mean that the operator acts on a covering manifold and is invariant 
under the unitary operators induced by the deck-transformations. 

Particular attention was devoted to the analysis of the gap structure of the spec- 
trum of periodic operators of Schrodinger type. More precisely, one is interested 
whether the spectrum in interrupted by spectral gaps, i.e. intervals on the real line 
which belong to the resolvent set. In case there are gaps: can one establish upper 
and lower bounds for the width and number of gaps and the spectral bands separat- 
ing them? Although the gap structure of the spectrum is a mathematically intrigu- 
ing question for its own sake, it is also important from the physical point of view. 
The features of gaps in the energy spectrum are relevant for the conductance prop- 
erties of the physical system. Even for periodic Schrodinger operators in Euclidean 
space the gap structure is highly non-trivial. This is maybe best illustrated by 
works devoted to the Bethe-Sommerfeld conjecture, e.g. |5591 12551 12551 I5571IT35] . 
Another interesting feature of some periodic Laplace-Beltrami operators is the ex- 
istence of L 2 -eigenfunctions, cf. the discussion in Remark 13.1.31 

These periodic operators on manifolds are generalised by their random ana- 
logues studied in this section. 

2.2. Random Schrodinger operators on manifolds: definitions. Let us 

explain the geometric setting in which we are working precisely: let A be a complete 
e?-dimcnsional Riemannian manifold with metric go- We denote the volume form 
of go by volo- Let T be a discrete, finitely generated subgroup of the isometries of 
(A, go) which acts freely and properly discontinuously on X such that the quotient 
M := A/r is a compact (<i-dimensional) Riemannian manifold. Let (0,£>n,P) be a 
probability space on which T acts by measure preserving transformations. Assume 
moreover that the action of T on f2 is ergodic. Now we are in the position to define 
what we mean by a random metric and consequently a random Laplace-Beltrami 
operator. 

Definition 2.2.1. Let {<7 w } w en be a family of Riemannian metrics on A. 
Denote the corresponding volume forms by voLj. We call the family {g^^en a 
random metric on (A, go) if the following five properties are satisfied: 
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(2.3) The map f2 x TX — > K, i— * guj(v,v) is jointly measurable. 

(2.4) There is a C g € ]0, oo[ such that 

C^ 1 3o(w, v) < g u (v, v) < C g g (v, v) for all v £ TX. 

(2.5) There is a C p <E]0,oo[ such that 

|V p w (a;)|o < C p for all a; G X, 

where Vo denotes the gradient with respect to go, p u is the unique smooth 
density of volo with respect to voLj, and \v\q = go(v,v). 

(2.6) There is a uniform lower bound [d — 1)K £ K for the Ricci curvatures of 
all Ricmannian manifolds (X, g u ). Explicitly, Ric(g a) ) > (d — l)Kg u for 
all ui £ f2 and on the whole of X. 

(2.7) The metrics are compatible in the sense that the deck transformations 

7: (X,g u ) ->• (X,g yu ), -f.x^~fx 
are isometries. 

Property 12.71) implies in particular that the induced maps 

U {u „y. L 2 (X,vol 7 - lw ) -» L 2 (X,voL), (tf (u , 7) /)(a;) = /(7 -1 *) 

are unitary operators. 

The density p^ appearing in l|2.5|l satisfies by definition 

/ f(x)dvol (x) = / f(x)p u (x)dvol u (x). 
Jx Jx 

It is a smooth function and can be written as 

Pa, (a;) = (detyo(eL,ei)) 1/2 = (det 5w (e , eg)) 

Here ej, . . . , e d denotes any basis of T X X which is orthonormal with respect to the 
scalar product go(x), and e* , . . . , e d € T^JT is any basis orthonormal with respect 
to g u (x). It follows from (|2.4|) that 

(2.8) Cg d/2 < Pu) (x) < C d g /2 for aU x £ X, u £ 
which in turn, together with property l|2.5l) and the chain rule, implies 

(2.9) |V Pot 1/2 (aOlo < /4 |V Pa,(a)|o for all x £ X, u £ fi 
Moreover, for any measurable A C X by 1)2. 8|) we have the volume estimate 

(2.10) C- d l 2 volo (A) < voL (A) < C d J 2 volo (A) 

We denote the Laplace-Beltrami operator with respect to the metric g w by A w . 

Associated to the random metric just described we define a random family of 
operators. 

Definition 2.2.2. Let {g w } be a random metric on (X, g ). Let V : £1 x X — > K 
be a jointly measurable mapping such that for all w € f2 the potential := 
■) > is in L\ oc (X). For each w G O let H u = — A w + K, be a Schrodinger 
operator defined on a dense subspace T> u of the Hilbert space L 2 (X, vol w ). The 
family {-ffo,}o,en is called a random Schrodinger operator if it satisfies for all 7 G T 
and u £ f2 the following equivariance condition 

( 2 - n ) = u (u,j)H^-i u U^ n ^ 
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Remark 2.2.3 (Restrictions, quadratic forms and selfadjointness). Some re- 
marks are in order why the sum of the Laplace-Beltrami operator and the potential 
is selfadjoint. We consider the two cases of an operator on the whole manifold X 
and on a proper open subset of X simultaneously. The set of all smooth functions 
with compact support in an open set A C X is denoted by (A) . For each u> G f2 
we define the quadratic form 

(2.12) Q(ff A ) : C C °°(A) x C C °°(A) -> R, 

(f,h) i * J^{x){vf{x),Vh{x)) dvol u (x) + J f(x)V u (x)h(x) dvo\ u (x) 

We infer from Theorem 1.8.1 in |68| that this quadratic from is closable and its 
closure Q(H%) gives rise to a densely defined, non- negative selfadjoint operator H^. 
Actually, Q(H^) is the form sum of the quadratic forms of the negative Laplacian 
and the potential. The result in |68| is stated for the Euclidean case X = W 1 but 
the proof works equally well for general Riemannian manifolds. 

The unique selfadjoint operator associated to the above quadratic form is called 
Schrodinger operator with Dirichlet boundary conditions. It is the Friedrichs exten- 
sion of the restriction H^\^joc /^y If the potential term is absent we call it negative 

Dirichlet Laplacian. 

There are special subsets of the manifold which will play a prominent role later: 

Definition 2.2.4. A subset T c X is called T -fundamental domain if it con- 
tains exactly one element of each orbit 0(x) := {y € X\ 3^ 6 Y : y = jx}, x € X. 

In j3J Section 3] it is explained how to obtain a connected, polyhedral T- 
fundamental domain T C X by lifting simplices of a triangularisation of M in a 
suitable manner. T consists of finitely many smooth images of simplices which can 
overlap only at their boundaries. In particular, it has piecewise smooth boundary. 

To illustrate the above definitions we will look at some examples. Firstly, we 
consider covering manifolds with abelian deck-transformation group. 

Example 2.2.5 (Abelian covering manifolds). Consider a covering manifold 
(X, go) with a finitely generated, abelian subgroup T of the isometries of X. If the 
number of generators of the group Y equals r, it is isomorphic to Z r ° x x . . . ZJ£ . 
Here = r and Z p is the cyclic group of order p. Assume as above that the 

quotient X/Y is compact. Periodic Laplace-Beltrami and Schrodinger operators on 
such spaces have been analised e.g. in 302, 2441 1245] , 

In the following we will discus some examples studied by Post in |244L 1245) . 
The aim of this papers was to construct covering manifolds, such that the corre- 
sponding Laplace operator has open spectral gaps. More precisely, for any given 
natural number N, manifolds are constructed with at least N spectral gaps. For 
technical reasons the study is restricted to abelian coverings. In this case the Flo- 
quet decomposition of the periodic operator can be used effectively. Post studies 
two classes of examples with spectral gaps. In the first case a conformal perturba- 
tion of a given covering manifold is used to open up gaps in the energy spectrum 
of the Laplacian. The second type of examples in |245) is of more interest to us. 
There, one starts with infinitely many translated copies of a compact manifold and 
joins them by cylinders to form a periodic network of 'pipes'. By shrinking the 
radius of the connecting cylinders, more and more gaps emerge in the spectrum. 
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Such manifolds have in particular a non-trivial fundamental group and are thus 
topologically not equivalent to R d . On the other hand their deck-transformation 
group is rather easy to understand, since it is abelian. In particular, it is amenable 
(cf. Definition 12 . 3 . 4|) . which is a crucial condition in the study conducted later in 
this section. Some of the examples in |244l I245| are manifolds which can be em- 
bedded in M 3 as surfaces. They can be thought of as periodic quantum waveguides 
and networks. See 244 for some very illustrative figures. 

Furthermore, in 245 perturbations techniques for Laplace operators on cov- 
ering manifolds have been developed, respectively carried over from earlier versions 
suited for compact manifolds, cf. |51l 1141 1114| . They include conformal perturba- 
tions and local geometric deformations. Floquet decomposition is used to reduce 
the problem to an operator on a fundamental domain with quasi-periodic boundary 
conditions and discrete spectrum. Thereafter the min-max principle is applied to 
geometric perturbations of the Laplacian. 

Related random perturbations of Laplacians are studied in |200l I198j , cf . also 
Example 12. 2. 71 In particular a Wegner estimate for such operators is derived. 

Now we give an instance of a covering manifold X with non-abelian deck- 
transformation group r. 

Example 2.2.6 (Heisenberg group). The Heisenberg group if 3 is the manifold 
of 3 x 3-matrices given by 



equipped with a left-invariant metric. The Lie-group if 3 is diffeomorphic to M 3 . Its 
group structure is not abelian, but nilpotent. 

The subset T = if3 n M(3, Z) forms a discrete subgroup. It acts from the left 
on if 3 by isometries and the quotient manifold H3/T is compact. 

Next we give examples of a random potential and a random metric which 
give rise to a random Schrodinger operator as in Definition 12.2.21 Both have an 
underlying structure which resembles alloy- type models (in Euclidean space). 

Example 2.2.7. (a) Consider the case where the metric is fixed, i.e. g u = go 
for all lu G Q, and only the potential depends on the randomness in the following 
way: 



Here u : X — > R is a bounded, compactly supported measurable function and 
q~f : Q — > K is a sequence of independent, identically distributed random variables. 
By considerations as in Remark 1 1 . 2 . 21 the random operator H u := — A + V U) u) G fi 
is seen to satisfy the equivariance condition. 

(b) Now we consider the situation where the metric has an alloy like structure. 
Let (go,X) be a Riemannian covering manifold and let a family of metrics {<? w } w 
be given by 



(2.13) 




(2.14) 




■yer 
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where u e C^°(X) and the r 7 : fl — >]0, oo[,7 E T are a collection of independent, 
identically distributed random variables. Similarly as in the previous example one 
sees that the operators A w are equi variant. 

2.3. Non-randomness of spectra and existence of the IDS. Here we 
state the main theorems on the non-randomness of the spectral components and 
the existence and the non-randomness of the IDS. They refer to random Schrodinger 
operators as defined in 12.2.21 

Theorem 2.3.1. There exists a subset f2' of full measure in (£!,,Efo,P) and 
subsets of the real line £ and £«, where • € {disc, ess, ac, sc,pp} such that for all 

a(H w ) = £ and a,(H u ) = S. 

for any • = disc, ess, ac, sc,pp. IfT is infinite, Y,disc = 0- 

The theorem is proven in |201) . see Theorem 5.1. The arguments go to a large 
part along the lines of 240J 11961 IT02 . Compare also the literature on almost 
periodic Schrodinger operators, for instance 298, 12lj. 

For the proof of the theorem one has to find random variables which encode the 
spectrum of {H^}^ and which are invariant under the action of T. By ergodicity 
they will be constant almost surely. The natural random variables to use are spec- 
tral projections, more precisely, their traces. However, since K is uncountable and 
one has to deal also with the different spectral components, some care is needed. 

Random operators introduced in Definition 12.2.21 are naturally affiliated to a 
von Neumann algebra of operators which we specify in 

Definition 2.3.2. A family {-B^} we n of bounded operators B u : L 2 (X, vo\ u ) — > 
L 2 (X, vol w ) is called a bounded random operator if it satisfies: 

(i) u> i— > (g^, Bufu;} is measurable for arbitrary f,g£L 2 (ilxX,¥o vol). 

(ii) There exists a w-uniform bound on the norms ||-Bo;|| for almost all u> € f2. 

(iii) For all u> G f2, 7 <E T the equivariance condition 

holds. 

By the results of the next paragraph § 12.41 {F(H U )} UJ is a bounded random 
operator for any measurable, bounded function F. 

It turns out that (equivalence classes of) bounded random operators form a 
von Neumann algebra. More precisely, consider two bounded random operators 
and {B^j}^ as equivalent if they differ only on a subset of £1 of measure 
zero. Each equivalence class gives rise to a bounded operator on L 2 (f2 x X, P o vol) 
by (Bf)(uj,x) := B LU f iAj (x), see Appendix A in |2(Jlj . This set of operators is a von 
Neumann algebra Af by Theorem 3.1 in |2(J1| . On J\f a trace r of type 11^ is given 

by 

t(B) := E [Tr( X rB.)} 
Here Tr := Tr w denotes the trace on the Hilbert space L 2 (X, voLj). Actually for 
any choice of u: f2 x X — > M + with X^er u -y- 1 Lj(l^ lx ) = 1 f° r au (^i^) E x X 
we have 



(2.15) 



t(B) = E [Ti(u.B.)} 
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In analogy with the case of operators which are T-invariant |19j we call r the T- 
trace. The spectral projections {P^ (] — oo, A[)} w of {H^}^ onto the interval ] — oo, A[ 
form a bounded random operator. Thus it corresponds to an element of Af which 
we denote by P(] — oo, A[). Consider the normalised T-trace of P 

(2 - 16) ^ (A) : = E[L.(^] 

The following is Theorem 3 in |199| . see also |201j . 

Theorem 2.3.3. PQ — oo, A[) is the spectral projection of the direct integral 
operator 



f 93 

Jn 



in 

and Nh is the distribution function of its spectral measure. In particular, the almost 
sure spectrum £ of {H u } u coincides with the points of increase 

{A € R| N H (X + e) > N H {\ - e) for all e > 0} 

ofN H . 

That the IDS can be expressed in terms of a trace on a von Neumann Algebra 
was known long ago. In 29?] and [298 Subin establishes this relation for almost 
periodic elliptic diffential operators in Euclidean space. He attributes the idea of 
such an interpretation to Berezin, see the last sentence in Section 3 of P 



We want to describe the self-averaging IDS by an exhaustion of the whole 
manifold X along a sequence A/ — > X, I € N of subsets of A. To ensure the existence 
of a sequence of subsets which is appropriate for the exhaustion procedure, we have 
to impose additional conditions on the group T. 

Definition 2.3.4. A group T is called amenable if it has an left invariant mean 

m L . 

Amenability enters as a key notion in Definition 12 .3 . 61 and Theorem 12. 3. 81 For 
readers acquainted only with Euclidean geometry, its role is motivated in Remark 

Under some conditions on the group amenability can be expressed in other 
ways. A locally compact group T is amenable if for any e > and compact K C V 
there is a compact G C Y such that 

m L (GAKG) < em L {G) 

where mj, denotes the left invariant Haar measure, cf. Theorem 4.13 in 242 . This 
is a geometric description of amenability of V. If T is a discrete, finitely generated 
group we chose mj, to be the counting measure and write instead | • | . In this case 
r is amenable if and only if a Foslner sequence exists. 

Definition 2.3.5. Let T be a discrete, finitely generated group. 

(i) A sequence {Ii}i of finite, non-empty subsets of T is called a Fxlner sequence 
if for any finite K C T and e > 

l/jAif/,1 <e|J,| 

for all I large enough. 
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(ii) We say that a sequence // C r, I £ N of finite sets has the Tempelman or 
doubling property if it obeys 

sup < oo 

(iii) We say that a sequence /; C T, / 6 N of finite sets has the Shulman property 
if it obeys 

IM-il „ 

sup < oo 

/eN Mil 

(iv) A Foelner sequence {Ii}i is called a tempered Foslner sequence if it has the 
Shulman property. 

In our setting T is discrete and finitely generated. (Actually, K :— {7 £ 
r I 7JF n T 7^ 0} is a finite generator set for T. This follows from the fact that 
the quotient manifold X/T is compact, cf. § 3 in [2].) Under this circumstances 
a Fcelner sequence exists if and only if there is a sequence J; C T, I £ N such 
that lim/^oo ^n^r^ = for all 7 € T. Moreover, for discrete, finitely generated, 
amenable groups there exists a F0lner sequence which is increasing and exhausts 
T, cf. Theorem 4 in [l]. 

Both properties (ii) and (iii) control the growth of the group T. Lindenstrauss 
observed in |21U| that each Fcelner sequence has a tempered subsequence. Note 
that this implies that every amenable group contains a tempered Foelner sequence. 
One of the deep results of Lindenstrauss' paper is, that this condition is actually 
sufficient for a pointwise ergodic theorem, cf. Theorem 12. 6. II Earlier it was known 
that such theorems can be established under the more restrictive Tempelman prop- 
erty |3()9l 11941 I31()| . Shulman 2^] first realised the usefulness of the relaxed 
condition (iii). 

In the class of countably generated, discrete groups there are several properties 
which ensure amenability. Abelian groups are amenable. More generally, all solv- 
able groups and groups of subexponential growth, in particular nilpotent groups, 
are amenable. This includes the (discrete) Heisenberg group considered in Example 
12.2.61 Subgroups and quotient groups of amenable groups are amenable. On the 
other hand, the free group with two generators is not amenable. 

For the discussion of combinatorial properties of Foelner sequences in discrete 
amenable groups see [J. 

Any finite subset I CT defines a corresponding set 

m :=int( (J 7^) 

where int(-) stands for the open interior of a set. 

In the following we will need some notation for the thickened boundary. Denote 
by do the distance function on X associated to the Riemannian metric g . For h > 0, 
let 9ftA := {x £ X\ da(%, <9A) < h} be the boundary tube of width h and A^ be the 
interior of the set A \ 9/jA. 

Definition 2.3.6. (a) A sequence {A/}; of subsets of X is called admissi- 
ble exhaustion if there exists an increasing, tempered Fcelner sequence {Ii}i with 
Ui Ii = r such that A z = fiilf 1 ), I £ N. 



INTEGRATED DENSITY OF STATES AND WEGNER ESTIMATES 



115 



(b) A sequence Ai,l G N of subsets of (X, go) is said to satisfy the Van Hove 
property 313] if 

(2.17) lim = o for all h > 
;^oo vol (A ; ) 

Remark 2.3.7. In our setting it is always possible to chose the sequences 
and {A;}; in such a way that they exhausts the group, respectively the manifold. 
However, this is not really necessary for our results. 

A simple instance where U;A; ^ X can be given in one space dimension. Let 
X = R, Z, J; = {0, . . . , I — 1}, J 7 = [0, 1] and consequently A/ = [0, 1]. One can use 
this sequence of sets to define the IDS of random Schrodinger operators although 
U/A; = [0, oo[. A non-trivial example where the sets Aj do not exhaust X can be 
found in 304, 306]. There Sznitman considers random Schrodinger operators in 
hyperbolic spaces, in which setting the approach presented here does not work due 
to lack of amenability. However, it is interesting that Sznitman obtains the IDS by 
choosing a sequence A/ which converges to a horosphere which is properly contained 
in the hyperbolic space. 

Thus an admissible exhaustion always exists in our setting. By Lemma 2.4 in 
243 every admissible exhaustion satisfies the van Hove property. Inequality (|2.1U|) 
implies that for a sequence with the van Hove property 

lim yoU^A) = Q for ^ h > Q 

holds for all u> G fi. Let us remark that one could require for the sets A; in the 
exhaustion sequence to have smooth boundary, see Definition 2.1 in |243j . Such 
sequences also exist for any X with amenable deck-transformation group L. This 
may be of interest, if one wants to study Laplacians with Neumann boundary 
conditions. For groups of polynomial growth it is possible to construct analoga 
of admissible exhaustions by taking metric open balls B ri (o) around a fixed point 
o G X with increasing radii n, . . . , r n , • • • — ► oo, cf. Theorem 1.5 in |243| . 

We denote by H l u the Dirichlet restriction of H u to A;, cf. Remark |2. 2. 31 and 
define the finite volume IDS by the formula 

Nl(X) := voUA;)- 1 ^ | X n (H l J < X} 

Now we are able to state the result on the existence of a self-averaging IDS. 

Theorem 2.3.8. Let {H^}^ be a random Schrodinger operator and T an ame- 
nable group. For any admissible exhaustion {A;}/ there exists a set Q' C f2 of full 
measure such that 

(2.18) lim fl£(A) = N H (X), 

l — >oo 

for every lu G Q' and every continuity point X G R of Nh ■ 

Definition 2.3.9. The limit in (|2.18|l is called integrated density of states. 



Thus all properties (i)-(iii) on page 1 1061 can be established for the model under 
study. In particular, formula (|2.18(l is a variant of the Pastur-Subin trace formula 
in the context of manifolds. The theorem is proven in §§ I2.4H2.6I It recovers in 
particular the result of Adachi and Sunada [2j on the existence of the IDS of periodic 
Schrodinger operators on manifolds. 
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Remark 2.3.10. Let us motivate for readers acquainted only with Euclidean 
space why it is natural that the amenability requirement enters in the theorem. 
In the theory of random operators and in statistical mechanics one often considers 
a sequence of sets A; , I G N which tends to the whole space. Even in Euclidean 
geometry it is known that the exhaustion sequence A;, I G N needs to tend to R d in 
an appropriate way, e.g. in the sense of Van Hove or Fisher 265 . Convergence in 
the sense of Van Hove |313) means that 

(2.19) lim J^M = 

i^oo I A; I 

for all positive e. 

If one chooses the sequence Ai,l € N, badly, one cannot expect the convergence 
of the finite volume IDS' N l u , I e N, to a limit. In a non-amenable geometry, any 
exhaustion sequence is bad, since (|2.19|l cannot be satisfied, cf. Proposition 1.1 in 

m 

Remark 2.3.11. We have assumed the potentials K, to be nonnegative and 
some of our proofs will rely on this fact. 

However, the statements of Theorem 12 . 3 . 1 1 on the non-randomness of the spec- 
trum and Theorem 12.3.81 on the existence of the IDS carry over to V u which are 
uniformly bounded below by a constant C not depending on u> G VI. Indeed, in 
this case our results directly apply to the shifted operator family {H^ — C} u en- 
This implies immediately the same statements for the original operators, since the 
spectral properties we are considering are invariant under shifts of the spectrum. 

2.4. Measurability. Since we want to study the operators as random 
variables we need a notion of measurability. To this aim, we extend the definition 
introduced by Kirsch and Martinclli 162 for random operators on a fixed Hilbcrt 
space to families of operators where the spaces and domains of definition vary with 

oj G n. 

To distinguish between the scalar products of the different L 2 -spaces we denote 
by (•, ')o the scalar product on L 2 (A, volo) and by || • ||o the corresponding norm. 
Similarly, and || • || w are the scaler product and the norm, respectively, of 

L 2 (A,voL). 

Definition 2.4.1. Consider a family of selfadjoint operators {H^}^, where the 
domain of is a dense subspace T> w of L 2 (A, voLj). The family {H^,}^ is called 
a measurable family of operators if 

(2.20) u^{f u ,F(H u )f u ) u 

is measurable for all measurable and bounded F : M — * C and all measurable func- 
tions / : £1 x A — > R with f(uj, ■) = f u G £ 2 (A, voLj) for every u € fi. 

Theorem 2.4.2. A random Schrodinger operator {i? w } w <=si as in Definition 
\2.2.£\ is a measurable family of operators. The same applies to the Dirichlet restric- 
tions {H^uen to any open subset A of X. 

For the proof of this theorem we need some preliminary considerations. 

Assumption l|2.4() in our setting implies that it is sufficient to show the weak 
measurability H2.2t)|l for functions / which are constant in u>. Note that I? (A, volo) 
and L 2 (A,voLj) coincide as sets for all ui £ fl, though not in their scalar products. 
Thus it makes sense to speak about / w = / G L 2 (A, voLj) " = " L 2 (A, volo). 
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Lemma 2.4.3. A random Schrddinger operator {Hu}u is measurable if and only 

if 

(2.21) lu (/, F(Hj)f) u is measurable 

for all measurable and bounded F: R — > C and all f G L 2 (A, volo). 

Proof. To see this, note that l|2.21ll implies the same statement if we replace 
f(x) by h(u>,x) — g{u>)f(x) where g £ L 2 (Vl) and / <E L 2 (A,volo). Such functions 
form a total set in L 2 (fl x A, P o vol). 

Now, consider a measurable h : £1 x A — > R such that h u :— h(u>, •) € L 2 (A, vol^) 
for every ui € il. Then h n (u>,x) :— Xh,n(oj) h(u>, x) is in L 2 (il x A,Povol) where 
Xh.n denotes the characteristic function of the set {u\ H/IwIIl^a.voI^) < n} C CI. 
Since Xh,n ~^ 1 pointwise on Q for n — > oo we obtain 

(hZ,F(H u )h») u - (h U) F(H u )h u ) u 

which shows that {i?^}^ is a measurable family of operators. □ 

To prove Theorem 12.4.21 we will pull all operators onto the same Hilbert 
space using the unitary transformation S u induced by the density p u 

S^. L 2 (A,vol ) -> L 2 (A,voL), (S u f){x) = p 1 J 2 (x)f(x) 

The transformed operators are 

(2.22) Au-^-S-^Su 

A,,; S- l V{^) C i 2 (A,vol ) — > i 2 (A,vol ) 

The domain of definition S'j 1 D(A^j) is dense in L 2 (A, volo) and contains all smooth 
functions of compact support in A. 

The first fact we infer for the operators A u , ui € f2 is that they are uniformly 
bounded with respect to each other, at least in the sense of quadratic forms. This 
is the content of Proposition 3.4 in |199j which we quote without proof. 

Proposition 2.4.4. Let Qq,Qu> be the quadratic forms associated to the oper- 
ators — Aq and A^, and T> c L 2 (A,vol ) the closure o/C^°(A) with respect to the 

norm (Qo(fJ) + \\f\\ 2 ) 1/2 - Then 

V = V(Q ) = V{QJ) 

and there exists a constant Ca such that 

(2.23) (Qo(fJ) + \\f\\l) < QM, /) + ll/llo < C A (Qo(fJ) + \\f\\ 2 ) ■ 
for all f S T> and uj E Q. 

Here T>(Q) denotes the domain of definition of the quadratic form Q. In the 
proof of this proposition the bound (|2.5|) — more precisely Q2.9[l — on the gradient 
of the density p^ is needed. It seems to be a technical assumption and in fact 
dispensable by using a trick from |69j . at least if A is precompact or of finite 
volume. 

Since we are dealing now with a family of operators on a fixed Hilbert space, 
we are in the position to apply the theory developed in |162j . The following result 
is an extension of Proposition 3 there. It suits our purposes and shows that our 
notion of measurability is compatible with the one in 162 . 
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Let Ti be a Hilbert space, T> C H. a (fixed) dense subset and £? w : I? — > 7i , G fi 
nonnegative operators. Denote by E = (J w a(B u ) the closure of all spectra, and by 
E c its complement. To establish the measurability of the family {B^}^ one can use 
one of the following classes of test functions: 

• Fx = {X]-oo,A[| A > 0}, 

• T 2 = {x ^ e itx \t E R}, 

• ^3 = ^^ e-^l t > 0}, 

• T A = {x i-> (z-x)- 1 ^ G C\E}, 

• J-5 = .7-4(2:0) = (zo — f° r a fixed zq E C \ E, 

• JF 6 = Ch = {/ : E — > C| / bounded, continuous}, 

• JF 7 = {/ : R — > C| / bounded, measurable}. 

The following proposition says, that it does not matter which of the above sets 
of functions one chooses for testing the measurability of {B^}^. 

Proposition 2.4.5. For i = 1, . . . , 7 the following statements are equivalent: 

(Fj) uj i ► (/, F(B UJ )h) -ft is measurable for all f,h € H and F G JFj, 

Proof. It is obvious that (F 4 ) => (F 5 ), (F 7 ) (F 6 ), and (F 6 ) (F 3 ). The 
equivalence of (Fi), (F2) and (F4) can be found in |162| . 
To show (F5) => (F4), consider the set 

Z := {z G E c | uj ^ (z — H^)~ l is weakly measurable } 

in the topological space E c . It is closed, since z n — > z implies the convergence of the 
resolvents, see e.g. |258l Theorem VI. 5]. A similar argument using the resolvent 
equation and a Neumann series expansion shows that z G Z implies B$(z) C Z 
where S := d(z, E). Since E c is connected, Z = E c follows. 

(F 3 ) (Fi): By the Stone- Weierstrass Theorem, see e.g. 258 Thm. IV.9], 
applied to C([0, 00]) it follows that JF 3 is dense in the set of Functions {/ G 
C([0, 00]) I /(oo) = 0} = Coo([0, oo[). We may approximate any X]_oo,a[ point- 
wise by a monotone increasing sequence < f n , n G N in Coo(E). Polarisation, the 
spectral theorem, and the monotone convergence theorem for integrals imply that 
X]-oo,\[{H u ) is weakly measurable. An analogous argument shows (Fi) => (F7), 
since any non-negative / G J-7 can be approximated monotonously pointwise by 
non-negative step functions /„, n G N. □ 

We use the following proposition taken from 295 (Prop. 1.2.6.) to show that 
{A^}uj is a measurable family of operators. 

Proposition 2.4.6. Let B^^uj G fi and Bq be nonnegative operators on a 
Hilbert space Ti.. Let Q^, uj G fi and Qq be the associated closed quadratic forms 
with the following properties: 

(2.24) Qu, uj G fi and Qo are defined on the same dense subset T> C H. 

(2.25) There is a constant C > such that 

C' 1 {Q (f, f) + \\f\\l) < QM, I) + ll/llo < C (Q (/, /) + WfWl) 

for all uj G fi and f G T>. 

(2.26) For every f E T> the map uj 1— ► Qu{f, f) is measurable. 

Then the family {B^}^ of operators satisfies the equivalent properties of Proposition 

rrot 
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By property (F7), this implies that {B^,}^ is a measurable family of operators. 

We apply the proposition to B^ = A u , where {4^}^ is defined in Proposition 
12.4.41 To do so we check that the properties l|2.4.6|l - (|2.2tj[l are satisfied: Since 
C£°(A) is dense in T>(Q U ) for all u>, the closures of C^°(A) with respect to any of 
the equivalent norms in (|2.23() coincide, which shows assumption l|2.4.6|) . Property 
is just J22HJ), property (j2~2^|l is obvious for / G C~(A) and follows by 
approximation for all / G T>. 

PROOF of Theorem 12.4.21 We already know that the transformed 'kinetic' 
part A w , w £ O of the Hamiltonian is measurable. To deal with the singular 
potential we introduce the cut off 

V"(x) := min{n, V u (x)} for n G N and u) G Q 

The auxiliary potential V™ is bounded and in particular its domain of definition is 
the whole Hilbert space L 2 (A,volo). Thus the operator sum 

AZ~A U + V£, coen 

is well defined and |162l Prop. 4] implies that it forms a measurable family of 
operators. To recover the unbounded potential V u , we consider the semigroups 
u> 1— > exp(— tA™),t > which are weakly measurable. 

The quadratic forms of ^4™ converge monotonously to the form of A^ :— A w + 
V u . Now Theorems VIII. 3. 13a and IX. 2. 16 in |152| imply that the semigroups of 
A™ converge weakly towards the one of A^ for n — > 00. Thus exp(— tA™) is weakly 
measurable, which implies the measurability of the family A™. 

Finally, since is multiplication with the measurable function (x, u>) 1— > p ul (x), 
this implies the measurability of the family H u = S^A^ S'J 1 , u> G fl. □ 

For later use let us note that the trace of measurable operators is measurable. 
More precisely we will need the fact that the mappings 

(2.27) w^Tr(xAe""^) and u ^ Tr(e~ tH ") 

are measurable. Note that one can chose an orthonormal basis for L 2 (A,voL_,) 
which depends in a measurable way on u>, cf. for instance Lemma II. 2.1 in |75j . 
Thus (|2.27|) follows immediately from the Definition l2.4.1l of measurable operators. 

2.5. Bounds on the heat kernels uniform in u>. This paragraph is devoted 
to heat kernel estimates of the Schrodinger operators Bl u . It consists of four parts. 
Firstly we discuss existence of L 2 -kernels of e~ tHl ° , t > and derive rough upper 
bounds relying on results in 68;. Secondly, we infer Gaufiian off-diagonal decay 
estimates of the kernels using estimates derived in |204| . We then present an idea 
of H. Wcyl to derive the principle of not feeling the boundary, and finally we state 
a proposition which summarises the information on the heat kernel needed in the 
next section. 

We have to control the dependence on the metric and potential of all these 
estimates since those quantities vary with the random parameter u> G f2. 

As is non-negative the semigroup e~ fH ",t > consists of contractions. 
Moreover, the semigroup satisfies some nice properties formulated in the following 
definition which enable us to derive estimates on the corresponding heat kernel. 
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Definition 2.5.1. Let A c X be open and /i a er-finite Borel measure on A. 
Let A be a real, non-negative, sclfadjoint operator on the Hilbert space L 2 (A,n). 
The semigroup e~ tA ,t > is called positivity preserving if e~ tA f > for any 
< / G L 2 (A, /i) and t > 0. Furthermore, e~ tA , t > is called a Markov semigroup, 
if it is well defined on i°°(A, fi) and the two following properties hold 

(2.28) e~ tA : L 2 (A, fi) — > L 2 (A, /i) is positivity preserving for every t > 

(2.29) e~ tA : L°°(A, //) -> L°°(A, /i) is a contraction for every t > 

In this case A is called a Dirichlet form. 

A Markov semigroup e~ tA is called ultracontr active if 

(2.30) e~ tA : L 2 (A, fi) -> L°°(A, /i) is bounded for all i > 

The above l|2.28ll and (I2.29|) are called Beurling-Deny conditions |2H 121] • 

We infer from 68 the following facts: A Markov semigroup is a contraction on 
L P (A, /u) for all 1 < p < oo (and all t > 0). For all w € ^ the Schrodinger operator 

on L 2 (A, voLj) is a Dirichlet form, |68l Thm. 1.3.5]. There the proof is given 
for X — M. d , but it applies to manifolds, too. By Sobolev embedding estimates and 
the spectral theorem e tA " is ultracontractive. Thus by Lemma 2.1.2 in 6§] each 
e tA ^ has a kernel, which we denote by k^, such that for almost all x, y € A 

(2.31) < kt{t,x,y) < ||e tA "|| 1)0o =: C A (t) 

Here ||B||i,oo denotes the norm of B : L 1 — * L 00 . For A = AT we use the abbreviation 
h x — h 

To derive an analogous estimate to l|2.31|l for the full Schrodinger operator with 
potential we make use of the Feynman-Kac formula. Using the symbol ~E X for the 
expectation with respect to the Brownian motion bt starting in x £ X the formula 
rCcicls 

(e- i ^/)( a; )=E x (e-/a^^)^/(M) 

For a stochastically complete manifold X and bounded, continuous V u the formula 
is proven, for instance, in Theorem IX. 7A in |93| . It extends to general non- 
negative potentials which are in L\ oc using semigroup and integral convergence 
theorems similarly as in the proof of Theorem X.68 in |255j . Since we consider 
(geodesically) complete manifolds whose Ricci curvature is bounded below, they 
are all stochastically complete, cf. for instance |131j or Theorem 4.2.4 in [136 . 

Since the potential is non-negative, the Feynman-Kac formula implies for non- 
negative / € i : (A, vohj) 

< (e^/) (x) < (e^ f ) {x) < ct{t) \\f\W 

for almost every a; € A. Thus e~ tH ^ : L 1 (A,vol a) ) — + L°°(A,vol w ) has the same 
bound C A (t) as the semigroup where the potential is absent. This yields the point- 
wise estimate on the kernel kg of e~ tH ^> : 

(2.32) < kff u (t, x, y) < C A {t) for almost every x, y e X. 

In the following we derive sharper upper bounds on the kernels which imply 
their decay in the distance between the two space arguments x and y. Such es- 
timates have been proven by Li and Yau 204 for fundamental solutions of the 
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heat equation. One would naturally expect that the fundamental solution and the 
L 2 -heat kernel of the semigroup coincide under some regularity assumptions. This 
is actually the case as has been proven for instance in |77| for vanishing, and in 
|199j for smooth, non-negative potentials. The proof in the last cited source uses 
that H^j is a Dirichlet form. 

To formulate the results of Li and Yau I204| which we will be using, we denote 
by : IxI-> [0, oo[ the Riemannian distance function on X with respect to g^. 
Note that the following proposition concerns the heat kernel of the pure Laplacian. 

PROPOSITION 2.5.2. For every t > there exist constants C(t) > 0, a t > 
such that 

(2.33) k u (t,x,y) < C(t) exp ( - a t dl(x, y)) 

for all co € Q and x, y G X . 

Proof. For a fixed Schrodinger operator the estimate (with do replaced by 
d u ) is contained in Corollary 3.1 in |204j . The properties l|2.4l) . I|2.8|l and 

C~ l d {x,y) <d w (x, y) < C g d Q (x,y) 

ensure that the constants C{t) and at in (|2.33(l may be chosen uniformly in uj. 
Moreover, for measuring the distance between the points x and y we may always 
replace d u by do by increasing at- □ 

Let us collect various consequences of Proposition 12.5.21 which will be useful 
later on. 

(i) The pointwise kernel bound on the left hand side of l|2.32|l can be chosen 
uniformly in lo € fi. 

(ii) We stated Proposition ^. 5. 2l for the pure Laplacian, although Li and Yau treat 
the case of a Schrodinger operator with potential. The reason for this is that 
we want to avoid the regularity assumptions on the potential imposed in 204 . 

To recover from <|2.33[) the case where a (non-negative) potential is present 
we use again the Feynman-Kac formula. We need now a local version of the 
argument leading to H2.32J) . more precisely we consider e~ tHl ° as an operator 
form L 1 (B e (y)) to L°°{B £ (x)) for small e > 0. Thus we obtain 

< k Ha (t, x, y) < C(t) exp ( - a t d 2 Q {x, y)) 

(iii) The estimates derived so far immediately carry over to the case where the 
entire manifold is replaced by an open subset Acl. 

< k^(t,x,y) < k H ^{t,x,y) 

This is due to domain monotonicity, see for example |68l Thm. 2.1.6]. 

(iv) The Bishop volume comparison theorem controls the growth of the volume 
of balls with radius r, see for instance |34j . (501 Thm. III. 6] or |42l . It tells 
us that the lower bound 12. 6|) on the Ricci curvature is sufficient to bound 
the growth of the volume of balls as r increases. The volume of the ball can 
be estimated by the volume of a ball with the same radius in a space with 
constant curvature K. The latter volume grows at most exponentially in the 
radius. For our purposes it is necessary to have an w-uniform version of the 
volume growth estimate. Using Properties (|2.4|) . (|2.ti|) and (|2.8|l we obtain the 
uniform bound 

voL ({y\ d u {x, y) < r}) < d e c ' 2r for all x G X 
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where Ci,C 2 do not depend on x and u>. This implies that for all exponents 
p > 0, there exists a M p (t) < oo such that the moment estimate 

f [k^Jt,x,y)r dvoUy) < M P (t) 

J A 

holds uniformly in A C X open, in a; € A and to € ft. We set M(t) := M\(t). 
(v) The heat kernel estimates imply a uniform bound on the traces of the semi- 
group localised in space. Let A C X be a (fixed) open set of finite volume. 
There exists a constant Cr r = Crr(A, t) > such that for all to G SI 

Tr( XA e- tff -) <C Tr 

Intuitively this is the same as saying that J A kn^ {t, x, x) dvol u (x) is uniformly 
bounded. However, since the diagonal {(a;, x)|a; € A} is a set of measure 
zero, the integral does not make sense as long as we consider fc# w as an L 2 - 
function. We do not want here to address the question of continuity of the 
kernel. Instead we use the semigroup property e ~ 2tHu ' = e~ "e~ tH ",t > 
and selfadjointness to express the trace as 
(2.34) 

Tr(x A e- tff ")= / f [k HLj (t/2,x,y)] 2 dvoUx)dvo\Uy) < Af 2 (i/2) voL(A) 

J A J A 

By (12.10(1 this is bounded uniformly in U! € ft. Applying domain monotonicity 
once more, we obtain 

(2.35) Tr(e-* H -) < M 2 (t/2) vol^A) 



The following lemma is a maximum principle for Schrodinger operators with 
non-negative potentials. Combined with the off-diagonal decay estimates in Propo- 
sition |^^21 it will give us a proof of the principle of not feeling the boundary. 

Lemma 2.5.3 (Maximum principle for heat equation with nonnegative po- 
tential). Let A C A fee open with compact closure, V be a non-negative func- 
tion, and u G C([0,T[xA) n C 2 (]0,T[xA) fee a solution of the heat equation 
JjU + (—A + V)u = on ]0, T[xA with nonnegative supremum s = sup{u(t, x) \ 
(t,x) e [0,T[xA}. Then, 

s = max < maxw(0, x), sup u(t, x) > 

[ x£A [0,T[x9A j 

Note that regularity of V is not assumed explicitly, but implicitly by the re- 
quirements on u. They are e.g. satisfied if V is smooth. Indeed, in that case the 
heat kernel is smooth, as can be seen following the proof of |68l Thm. 5.2.1]. 

Now we are in the position to state the second refined estimate on the heat 
kernels, the principle of not feeling the boundary. It is a formulation of the fact 
that the heat kernel of the Laplacian restricted to a large open set A should not 
differ much from the heat kernel associated to the Laplacian on the whole manifold, 
as long as one stays away from the boundary of A. As before, we derive this 
estimate first for the pure Laplacian and then show that it carries over to Schro- 
dinger operators with non-negative potential. 
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PROPOSITION 2.5.4. For any fixed t, e > 0, there exists an h = h(t, e) > such 
that for every open set Acl and all u> £ Q 

Q<k u {t,x,y)-kt{t,x,y) < e, 

for all x £ A, y £ A^. 

Proof. The first inequality is a consequence of domain monotonicity. So we 
just have to prove the second one. 

Fix well and t,e > 0. Choose h > such that 

v2> 



C(t) exp(-a t (V2) : 



< e 



Note that the choice is independent of u>. For any y £ A/j and < 5 < h/2 
denote by Bg(y) the open d -ball around y with radius S. Let fs € C^°(Bs(y)) be 
a non-negative approximation of the ^-distribution at y. 

We consider now the time evolution of the initial value / under the two semi- 
groups generated by A w and A A , respectively. 

U!(t,x):= / k u (t, x,z)f s (z)dv6l w (z) = / k u (t,x,z)fs(z)dval u> (z). 
Jx J A 

u 2 {t,x) := / k^(t,x,z)fs{z)dvol u (z). 



The difference x) — u^it, x) solves the heat equation J^-u = A^u and satisfies 
the initial condition Wi(0, x) — 1*2(0, x) — f$(x) — fg(x) — for all x £ A. Now, by 
domain monotonicity we know k^(t, x, z) — k^(t, x, z) > 0, thus 

ui(t,x) - u 2 (t,x) = / [k u (t, x,z) - fc A (i,x, z)] f$(z) dvolu(z) > 



for all t > and x £ A. The application of the maximum principle yields 

(2.36) u\(t,x) — U2(t,x) < max {ui(s, w) — u?,(s, w)\ . 

]o,t]x9A 

The right hand side can be further estimated by: 

u\(s,w) - v,2(s,w) < / k u} (s ) w,z)fg(z)dvo] w (z) = / k u} (s,w,z)fs(z)dvo] w (z). 

JA JA h/2 

Since w £ dA and z £ A h / 2 , we conclude using Proposition 12 . 5 . 'J\ 

k u (s,w,z)fs(z)dvo\ w (z) < C(t)exp(~ a t {h/2) 2 ) < e 



Since the bound is independent of 5 we may take the limit 8 — > which concludes 
the proof. □ 

One can prove the principle of not feeling the boundary by other means too, 
see for instance [2151 1791 12431 . This alternative approach uses information on 
the behaviour of solutions of the wave equation. Unlike the solutions of the heat 
equation, they do not have the unphysical property that their support spreads 
instantaneously to infinity. Actually, the solutions of the wave equation have finite 
propagation speed |307) . Fourier transforms and the spectral theorem turn this 
information into estimates on the difference of the solutions of the free and restricted 
heat equation. Sobolev estimates lead then to the principle of not feeling the 
boundary. 
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See also Section 7 in [254 . 

Remark 2.5.5. Similarly as in Lemma 12.5.31 one can prove the proposition 
if a potential is present. More precisely, Proposition 12.5.41 is valid for Schrodinger 
operators with potentials V such that for continuous initial and boundary values the 
solution of the heat equation j^u = -(-A u +V)u is in C([0,T[xA) n C 2 (]0,T[xA). 
However, Proposition 12.5.41 implies an analogous estimate for the case where a non- 
negative potential is present, similarly as in (0 on page ll2ll 

Consider e~ tH ^ — e~ tH -> as an operator from L 1 (A/i) to L°°(A), and denote 
by the first exit time from A for a Brownian motion starting in x. By the 
Feynman-Kac formula, we have for < / € L 1 (A? l ) 

[(e-**" - e~^)f]{ X ) = E x (e-/o W /(6t)x{6|r ^ t} ) 

< E x ^/(6 s )x{6| T A< t} ^ = [[k u (t,x,y)-k£(t,x,y)]f(y)dvo] w < e J f(y)dvo\ u 

if we chose h as in Proposition 12 . 5 . 41 Thus for almost all x S A, y £ Ah 
(2.37) k Hu (t,x,y)-k^Jt,x,y) < \\ e - tH " - e~ tH " \\ L i (Ah )^L°°(A) <e 



The upper bounds on the heat kernel and the principle of not feeling the bound- 
ary enable us to prove a result on the traces of localised heat-semigroups. 

In the macroscopic limit, as A tends (in a nice way) to the whole of A, the two 
quantities 

Tr(x A e-* H ") and Tr( e -* K *) 
are approximately the same. The precise statement is contained in the following 

Proposition 2.5.6. Let {A;}/ s n, be an sequence of subsets of X which satisfies 
the van Hove property \2. 1 1\ and let {H^}^ be a random Schrodinger operator. Then 



lim sup — - 

Z^oo^gQ VOL (A;) 



Tr( XAi e-^)-Tr( e -*"'. 



Q 



Proof. We consider first a fixed I € N and abbreviate A = A/. For the 
operator e~ tH ^ we may write the trace in the same way as in Ij2.34|l to obtain 

(2.38) Tr(e- tH *)= f [ [fc& (t/2, x, y)] 2 dyol u (x)dvol u (y) 

J A J A 

We express the difference of (|2.34fl and 12.38[) using 

(k H J 2 - (4 J 2 - (k Hul - kbj{k Hu + 4j 

Next we chose h = h(t/2,e) > as in Proposition 12.5.41 and decompose the inte- 
gration domain according to 

A x A = (A x Ah) U (A x d h A) 
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The difference of the traces can be now estimated as 

(2.39) Tr( X A e~ tH " ) - Tr(e~ tH » ) 

f f t t t t 

= JaJa [ kH "(~2 ,X,y ^ ~ kH ^2' X ' V ^ [ kH ">(~2 ,X,y ) +kH »(~2 ,X,y ^ dY0 ^ x ^^ 

nt t t t 

| a [kH u (^,x,y) -k Hai (-,x,y)] [k Hul (^,x,y) +k^ u (-,x,y)] dvo\,(x,y) 

The first term is bounded by 2M(t/2) e voL(A) and the second by 

2M{t/2)C(t/2)vol w (d h A) 

It follows that 

< ^A)( Tr (XA^^) -Tr(e^)) < 2M(t/2)e + 2M(t/2)C(t/2) V °|^) 

Now, we let I go to infinity. Since the sequence A; satisfies the van Hove property 
(|2.17|l and since our bounds are uniform in ui, the proposition is proven. □ 

2.6. Laplace transforms and Ergodic Theorem. This paragraph com- 
pletes the proof of Theorem l2 . 3 . 81 It relies, apart from the results results established 
in the previous paragraphs 12. 41 - 12. 51 on a general ergodic theorem and a criterion 
for the convergence of distribution functions. 

Lindenstrauss proved in |21Ul I2U9| an ergodic theorem which applies to locally 
compact, second countable amenable groups. It includes as a special case the 
following statement for discrete groups. 

Theorem 2.6.1. Let T be an amenable discrete group and (£l,Bn,¥) be a prob- 
ability space. Assume that T acts ergodically on Q by measure preserving transfor- 
mations. Let {Ii}i be a tempered Foslner sequence in T. Then for every f 6 L x (f2) 

(2.40) lim Ir £/(7w)= E(/) 

for almost allw£ f2. 

In the application we have in mind / € L°° , so the convergence holds in the 
L 1 -topology, too. 

Remark 2.6.2. Some background on previous results can be found for instance 
in Section 6.6 of Krengel' s book [194] . in Tempelman's works 308] 13091 1310] or 
some other sources |95H15ll94l l236 . The book 310 gives in § 5.6 a survey of Shul- 
man's results |276| . Mean ergodic theorems hold in more general circumstances, 
see for instance [1941 § 6.4] or [3101 Ch. 6]. 

We will apply the ergodic theorem above not to the normalised eigenvalue 
counting functions N l u , but to their Laplace transforms C 1 ^. The reason is, that the 
C l u are bounded, while the original Nfj are not. The following criterion of Pastur 
and Subin |237L 1298) says that it is actually sufficient to test the convergence of 
the Laplace transforms. 

Lemma 2.6.3 (Pastur-Subin convergence criterion). Let N n be a sequence of 
distribution functions such that 
(i) there exists a Aq € R such that Ni(X) = for all A < Aq and I G N, 



126 



I. VESELIC 



(ii) there exists a function C : K+ -> K such that C l (t) := J e- xt dNi(X) < C(t) 
for all I € N and < > 0, 

(iii) lirn;_ >(X) =: C(t) exists for all t > 0. 

TTien £ is the Laplace transform of a distribution function N and for all continuity 
points A of N we have 

N(X) := lim TV/ (A) 

l — >oo 

Finally, we present the proof of Theorem 12.3.81 on the existence of a self- 
averaging IDS: 

Proof of Theorem 12.3.81 We have to check the conditions in the previous 
lemma for the normalised eigenvalue counting functions . The first one is clearly 
satisfied for Ao = 0, since all operators we are dealing with are non-negative. To 
see (ii), express the Laplace transform by the trace of the heat semigroup 

uK ' voL A voL (A V ; 

The sum extends over all eigenvalues A„ of H 1 ^, counting multiplicities. Now, (|2.35() 
implies condition (ii) of the Pastur-Subin criterion. 

To prove (iii) we will show for all t > the convergence 

lim Cl(t) = f e- tx dN H (\) 

111 (i 1 and) P almost sure-sense. 

Introduce for two sequences of random variables ai(u>), bi(u>), I € N the equiva- 
lence relation ai 3 bi if they satisfy a; — bi — > 0, as I — > oo, in L 1 and P-almost 
surely. 

For technical reasons we will deal separately with the convergence of the enu- 
merator and denominator in 

£l ) W=voL(A / )- 1 Tr(e- t ^) 

However, we need some normalisation, to avoid divergences. Consider first the 
enumerator with an auxiliary normalisation 

(2.41) \I l \- l Tr(e-* H, ») 

By Proposition 12 . 5 .fil the equivariance, and Lindenstrauss' ergodic theorem l2.fi. II 

75/r 1 

= l/r 1 £ Tr(x^ e -* H -) 3 ~°° E{Tr( X ^e- tff -)} 

76-fi 

Similarly we infer for the normalised denominator 

\h\-\o\,{h) = £ voL( 7 ^) = \It\~ 1 £ vol 7W (^) j ~°° E {vol.(.F)} 
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Note that by (|2.1()|1 all terms in the above line are bounded form above and below 
uniformly in ui. By taking quotients we obtain 

, \h\-^T{e-^) ^oo E {Tr(xre- tH ')} 

u{> voL (A,) ~ E{vol.(F)} 

The right hand side is the Laplace transform of Nh, see the proof of Theorem 6.1 
of 201 for a detailed calculation. □ 

2.7. Approach using Dirichlet-Neumann bracketing. We sketch briefly 
an alternative proof of the existence of the IDS due to Kirsch and Martinelli |161| . 
It applies to random Schrodinger operators on M d . 

It relies on a ergodic theorem for superadditive processes by Akcoglu and Kren- 
gel |12| and estimates on the number of bound states essentially implied by the Weyl 
asymptotics. 

Let us explain the notion of an superadditive process in our context. Denote by 
Z the set of all multi-dimensional intervals or boxes A in M d such that A = {x\ cij < 
Xj < bj, for j = 1, ... , d} for some a, b G Z d with a,j < bj for all j = 1, . . . , d. The 
restriction of H u to an A G Z with Dirichlet boundary conditions is denoted by 

and with Neumann boundary conditions by H^' N . Consider a group {Tfc} fegZ d 
(or semigroup {7fc} fceN d) of measure preserving transformations on the probability 
space (fi,B n ,P). 

Definition 2.7.1. A set function F: Z — ► L 1 ^) is called a (discrete) super- 
additive process (with respect to {Tk}k) if the following conditions are satisfied 

(2.42) F A o T k = F A+fc for all k G Z d (or N„) ,AeZ 

n n 

(2.43) if Ai, . . . , A n G Z such that A := int ( (J A fc ^) G Z then F A > ^ F Afc 

fe=i fe=i 

(2.44) 7 := 7 (F) := sup | A^ 1 E {F A } < oo 

Aez 

F is called subadditive if — F is superadditive. 

Similarly one can define continuous superadditive processes with respect to an 
action of M. d on Q. 

We formulate the main result of |12| in the way it suits our needs. 

Theorem 2.7.2. LetF be a discrete superadditive process and hi := [— 1/2, l/2] d , 
I G N. Then the limit lim;—^ l~ d F\, exists for almost all uj G O. 

More generally, one can replace the A;, / e N by a so called regular sequence, 
cf. USEE H21 QUI or § 6.2 in [T31| . 

To apply the superadditive ergodic theorem we consider for arbitrary, fixed 
Asl the processes given by the eigenvalue counting functions of the Dirichlet and 
Neumann Laplacian 

Ff :=F°(\,u;)~#{n\\ n (H^<\}, KeZ 
Fa ■= ^A (A,w) := #W X n (H^ ,N ) < X}, AeZ 

where H u is a random operator as in Definition 11.2.31 Obviously for A = A/ = 
[-l/2,l/2] d we have Ff (A) = l d Nl{\). We will show that Ff ,A G Z is a super- 
additive process, which is also true for — F A , A G Z. Property (|2.42|l follows from 
the equivariance of {F^}^, while (|2.43|l and Ij2.44|l are implied by the following 
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Lemma 2.7.3. Let be a random operator as in Definition U.2.31 and A a 
fixed energy value. 

(i) For two cubes Ai C A 2 we have > Ffc and F^ > Ffc . 

(ii) If Ai, A 2 G Z are disjoint such that A = Ai U A2 U M G Z where M C R d is 
a set of measure zero, then 

F? > F£ + Fg 

f?<fZ+f£ 

(iii) There exists an C\ G M. such that for all A £ Z and uj £ Q we have F® (uj) < 
C A |A|. 

Proof. The first two statements are known as Dirichlet-Neumann bracketing 
and are stated e.g. in Proposition XIII. 15.4 in |256j . See also Section 1.5 in |50j 
for analogous results on manifolds. 

To prove (iii) we consider first the case where the potential is identically equal 
to zero and infer from the Weyl asymptotics that there is a constant C(A) G K such 
that 

#{n|A„(-A A )<A}<C(A)|A| 

Here A A denotes the Laplace operator on A with Dirichlet boundary conditions. 
Since the potentials we consider are infinitesimally bounded with respect to the 
Laplacian, uniformly in u, the Min-Max principle for eigenvalues implies the same 
estimate for the full Schrodinger operator, with an suitably changed constant. See 
Sections XIII. 3, 15 and 16 in 256 for more details and, for explicit estimates, the 
works |264L I206L I205L 164) deriving the Lieb-Thirring and Cwikel-Lieb-Rozenblum 
bounds on the number of bound states. □ 



Now we can state the main result of [161] . 

Theorem 2.7.4. There exists a set f2' C fl of full measure such that 
(2.45) N{\) := lim N l u (X) 

l — >oo 

exists for every u> G Of and every continuity point A G M of N . 

Proof. For a fixed A G R one applies the ergodic theorem of |12| to the 
superadditive process F D — F D (A,ui). Since in our case the transformation group 
is ergodic, the limit N(X) equals r y(F D (X)), in particular it is independent of u>. 
Almost sure convergence means that there is a set £l\ of measure one for which the 
convergence holds. Let S C M. be dense countable set. Then f2' = Dx^s^x still 
has full measure and (|2.45l) holds for all A G S and u G CI'. Since S was dense, 
this shows the convergence (|2.45|) at all continuity points of the limit function. 
Afterwards one modifies the limit function N to be left continuous. See |161l 115?) 
for details. □ 



For models which satisfy both the conditions of the previous theorem and of 
the two IDS's coincide. 

Under certain regularity conditions the theorem remains true if Neumann bound- 
ary conditions are used to define the IDS, and also for M d -ergodic potentials, cf. for 
instance [16T1H42] . 
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2.8. Independence of the choice of boundary conditions. Consider again 
the more general setting of Schrodinger operators on a Riemannian covering man- 
ifold X . If the open subset A C X of finite volume is sufficiently regular, the 
Neumann Laplacian H^' N on A has discrete spectrum. One condition which en- 
sures this is the extension property of the domain A, see e.g. |68| . which is in turn 
satisfied if the boundary <9A is piecewise smooth. Minimal conditions which ensure 
the extension property are discussed in § VI. 3 of |292| . Thus it is possible to define 
the normalised eigenvalue counting function 

N^ N (X) := i #{n e N\ K(H^ N ) < A} 

Let A; be an admissible exhaustion Aj C 1,1 £ N of sets which all have the 
extension property. Consider the sequence of distribution functions Njj N := N^ l ' N . 
It is natural to ask whether it converges almost surely, and, moreover, whether its 
limit coincides with N as defined in 1)2.45)1 . If this is true, the IDS is independent of 
the choice of Dirichlct or Neumann boundary conditions used for its construction. 
This indicates that boundary effects are negligible in the macroscopic limit. 

However, this turns out not to be true for all geometric situations. Sznitman 
studied in |304L 1306) the IDS of a random Schrodinger operator on hyperbolic 
space with potential generated by a Poissonian field. He showed that the IDS does 
depend on the choice of boundary condition used for its construction. Actually, 
he computes the Lifshitz asymptotics of the IDS at energies near the bottom of 
the spectrum and shows that it is different for Dirichlet and Neumann boundary 
conditions. 

In contrast, in the case of Euclidean geometry X — M. d , the question of bound- 
ary condition independence has been answered positively already some decades ago 
[271 11611 12991 188| for a large class of Z d or R d -ergodic random potentials. Re- 
cently, there has been interest in the same question if a magnetic field is included in 
the Hamiltonian, see also § 15.61 In this case the coincidence of the IDS defined by 
the use of Dirichlet and Neumann boundary conditions was established for bounded 
potentials in 232 , for non- negative potentials in |82) . and for certain potentials 
assuming both arbitrarily large positive and negative values in |142j and |139j . 
The last mentioned approach seems to be extensible to non-Euclidean geometries 
which is matter of current research. 



3. Wegner estimate 

In 1981 Wegner 325 proved on a physical level of rigour a lower and upper 
bound on the density of states (DOS) of the (discrete) Anderson model and similar 
lattice Hamiltonians. The density of states — if it exists — is the density function of 
the IDS. Wegner's argument did in particular not rely on any information about the 
type of the spectrum in the considered energy interval. This was important since 
before Wegner's result there where various conjectures in the physics community 
how the DOS should behave at the mobility edge, if it exists. This is the name for 
the critical energy which is supposed to form the boundary between an interval with 
pure point spectrum and another one with continuous spectrum. Note however that 
there is so far no rigorous proof of the existence of continuous spectrum for ergodic 
random Schrodinger operators. 
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There were arguments suggesting that the DOS should diverge to infinity at the 
mobility edge, others that it should vanish. Wegner's estimate discarded this mis- 
conceptions. In a sense it is a negative result: you cannot recognize the borderline 
of different spectral types by looking at the IDS. 

In the sequel we concentrate on alloy type models as defined in 11.2.11 (and 
Remark 11.2. 2|1 . We will be concerned here with upper bounds on the DOS only. 
It is derived by considering its analoga on finite boxes. So what we are speaking 
about in this section is an estimate on 

E {TrP'(J)} = |Af|E {N l u (E 2 ) - A^i)} 

where for the moment for notational convenience we only consider half open energy 
intervals / = [Pi , E 2 [. By the Cebysev-inequality one sees 

(3.1) P{<r(#i) n I £ 0} = P{TrP<(J) + 0} < E{TrP< (/)} 

This means that a bound on the averaged trace of the projection gives immediately 
a bound on the probability to find an eigenvalue in the considered energy interval. 
Actually, in the literature on Anderson localisation often the (weaker) bound on 
P {TrP^(I) 7^ 0} is called Wegner estimate, since it is sufficient for the purposes of 
multiscale analysis, see § 13.21 

In the following we will adopt the following notation: A; :=] — l/2,l/2[ d C M. d 
denotes the cube of side length I centred at zero. Occasionally we suppress the 
dependence on the size and just write A. A cube centred at x £ M. d is denoted by 
A; + x — {y + x\y £ A;} or Aj(x). The characteristic function of the unit cube 
Ai + j is abbreviated by Xj- For I £ N the symbol A; denotes the lattice points in 
Z d such that 

int ( y AT(7)) = A, 

More explicitly: A ; = A; n Z d . 

3.1. Continuity of the IDS. The estimates on the expected number of en- 
ergy levels in /, which most authors derive or use (for localisation proofs) are 
'polynomial', more precisely 

(3.2) E{TrP<(/)} <C w \I\ a \Ai\ b 

Here |/| and |A;| denote the (1-dimensional, respectively d-dimensional) Lebesgue 
measure of the energy interval J, or the set A;, respectively. The Wegner constant 
Cw depends on the various parameters of the model and for continuum Hamil- 
tonians on the supremum of I. Actually, Cw can be assumed to be a monotone 
non-decreasing function of | sup/|. However, once sup / is fixed, Cw is independent 
both of the energy interval length and the volume. It is obvious that the energy 
and volume exponents must satisfy a £]0, 1], b £ [1, oo[. As far as the exponents are 
concerned, the Wegner estimate is optimal if the dependence on the volume and 
energy length is linear, i.e. a = b = 1. 

For, if b = 1, the bound l|3.2|l carries over to the infinite volume IDS: 

(3.3) lim E {Nl{E 2 ) N l u (E 2 )} = km gi^ffiMl < C ' w \E 2 - E 1 \ a 

l—>oo I— >oo An 
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Since we know from Theorem 12.3.81 and dominated convergence that for E\ , E 2 in 
a dense set of energies 

lim E {Nl(E 2 ) - N l u (Ei ) } = N(E 2 ) - N(E 1 ) 

it follows 

N(E 2 ) ~ N(Ei) < C w \E 2 - E^ 

where Cw = Cw(E 2 ). Now the monotonicity of the IDS implies its Holder conti- 
nuity. Moreover, if the estimate l|3.2|) is linear in the energy, i.e. a = 1, the IDS is 
even Lipschitz continuous. Thus its derivative, the density of states, 

exists for almost every energy figl and is locally bounded 

n(E) < C W (E 2 ) for all E < E 2 

So the Wegner constant turns out to be a locally uniform bound on the DOS. 

Remark 3.1.1. For certain models the bounds derived on P {cr(iJ^) n / ^= 0} 
are not polynomial in the volume. This is the case for one-dimensional Anderson or 
alloy type models where the coupling constants uij , j £ 1 are distributed according 
to the Bernoulli distribution: for some p £E]0, 1[ the random variable loq takes on 
the value 1 with probability equal to p and the value with probability 1—p. Since 
this disorder regime is highly singular, the 'usual' proofs of the Wegner estimate, 
see Sections 0] and 13 fail. The ones that do work, yield somewhat weaker results. 
Namely, it is proven in j46| (cf. also 275 ) for the discrete Anderson model and in 
66 for the continuum alloy type model, that for a fixed compact energy interval 
/ and all (3 G]0, 1[, 7 > there exist lo £ N and a > such that 

(3.5) P{d(<j(H l J,E)<e-^}<e- al " 

for all E G / and all I >Iq. Here in the case of the continuum model it has to be 
assumed that / is disjoint from a discrete set of exceptional energies. 

This estimate obviously does not imply a continuity estimate for the infinite 
volume IDS. Interestingly, for these models, the Holder continuity of the IDS is 
established using other techniques which are specifically tailored for the one dimen- 
sional case, see |197| . |46l App. to § 5] and 66 Thm. 4.1]. Subsequently, this 
regularity result is used to derive the finite volume estimate (|3.5|) . In higher di- 
mensions, as we have discussed above, one proceeds in the other direction, carrying 
over finite volume estimates to the macroscopic limit. 

The bound (|3.5(l is still sufficient as an input for the multiscale analysis which 
proves localisation, cf. e.g. [52~0] or our discussion in § 13.21 In the discussion of 
Spencer's example in S 13.21 we will obtain an insight why subexponentially small 
eigenvalue splittings are effective enough to prevent resonances. 

Remark 3.1.2 (Continuity of the IDS on the lattice and in one dimension). In 
|74j Delyon and Souillard showed by a very simple argument that the IDS of the 
Anderson model on the lattice Z d is continuous, regardless of the continuity of the 
distribution of the potential values. The potential may even be correlated over long 
distances, as long as it is an ergodic stochastic field. Delyon and Souillard use the 
unique continuation property of the discrete Schrodinger equation, to prove that no 
eigenvalue can be sufficiently degenerated to produce a jump of the IDS. At the end 
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of Remark 13 . 1 . 31 we contrast their theorem with the situation in graphs other than 
the lattice. In |62| a quantitative version of the continuity is proven: for random, 
ergodic lattice Hamiltonians the IDS is actually log-Holder continuous. See also 
|212| and [1011 1801 1531 12201 181| for a related result for graph Hamiltonians. 

Similarly, the IDS is continuous for one dimensional Hamiltonians, both on Z 
and on R, |240L 11491 121j . Again, this result can be strengthened to log-Holder 
continuity, cf. |63| . 

Remark 3.1.3 (Continuity of the IDS and geometry). So far we have only 
mentioned proven or expected assertions on the continuity of the IDS. One might 
ask whether there are interesting models which exhibit a discontinuous IDS. It 
turns out that this phenomenon may occur, if the configuration space has a more 
complicated geometry than 1 d or W l . Another example would be the IDS of the 
Landau Hamiltonian, cf. e.g. the references in § 15.61 in particular 232 . 

Maybe the simplest example to illustrate the difference between Euclidean and 
more general geometry is provided by periodic Schrodinger operators. Under mild 
assumptions on the Z d -periodic potential V P er the IDS of the Schrodinger operator 
H = -A + Vp Cr on R d is absolutely continuous, cf. e.g. [2561 Problem 145]. |3"TT] . 
In particular the IDS cannot have jumps. However, precisely this can occur for 
Laplace-Beltrami operators (even without potential) on a Ricmannian covering 
manifold A, as it was mentioned in |301l App. 2] and is a subject of current research 
198, 200 . This phenomenon can be deduced from the fact that Laplacians on 
covering manifolds may have eigenvalues, as has been shown in 186 . Furthermore, 
the size of the jumps of the IDS is related to certain geometric invariants. Examples 
of such invariants are the order of the torsion subgroup of the deck transformation 
group r of X and the L 2 -Betti numbers of A, which can be expressed in terms of the 
L-trace on a certain von Neumann algebra, establishing the connection to the repre- 
sentation of the IDS discussed before TheoremES3| 751 1791150112^91181] . Related 
material can be found in QH 13021 12121 l300l 12681 11301 12701 l220l fTSl 12131 1214] . 

Some of the Wegner estimates we present in Sections 0] and extend to alloy 
type models on manifolds. A particularly interesting phenomenon occurs if one con- 
siders a periodic Laplace-Beltrami operator with discontinuous IDS, and perturbs it 
randomly such that the IDS of the perturbed operator is continuous. This happens 
if either an appropriate alloy type potential is added to the Hamiltonian or if the 
metric is multiplied by an appropriate alloy type perturbation, see 1981 1200] , 

A discontinuous IDS may also occur for models with a random geometry. This 
is the case for the tight-binding Hamiltonian defined on Delone sets studied in 
[27T21IT75] . Ideas related to the ones in [175] have been used in |2"B"7| and in § 2 
of |220| . The paper |267| is devoted to the proof of the existence of spectral gaps 
for certain graph Hamiltonians. 

We will discuss a different example, the quantum percolation model, in some 
more detail, since it fits readily in the class of models which we have described so far. 
This model has been studied amongst others in [721 [7IJ 11551 12721 1521 1314| . We 
sketch the site percolation problem on 1 d with probability parameter p above the 
percolation threshold: let Vk : fl — > {0, oo},£; £ 1 d be a sequence of independent, 
identically distributed random variables which take on the value with probability 
p and the value oo with probability 1 — p. Define A w to be the infinite component 
of the set of active sites {k e Z d \ Vk(oj) = 0}. The graph X u is called the (active) 
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infinite cluster. For p above a certain critical value p c it is known that almost surely 
an infinite cluster exists |154L Il32] and is unique [51 1116) . 

One defines the Laplacian h u on as the restriction of the finite difference 
operator onto ^(X^). For a cube A; one defines X l u to be those active sites in A/nZ d 
which are connected to the boundary <9A; by a chain of active sites. The finite 
volume Laplacian is the usual finite difference operator restricted to Z 2 (X^). 

Although the finite active clusters, which would obviously give rise to bound 
states, are not taken into consideration, it turns out that h u has bound states. This 
was as first observed in 155 . Eigenstates with finite support in the infinite cluster 
are called molecular states. The existence of such states affects the properties of 
the IDS of h^, which is defined in the following way. For each I £ N the normalised 
eigenvalue counting function of the Hamiltonian h 1 ^ is defined as 

KiE) ■■= W ^ W) #in\KiK)<E} 

which converges for I —* oo to an non-random limit almost surely 52J 1314) . A 
simple construction shows that there are locally supported eigenf unctions, which 
depend only on the pattern of X^ in a bounded region. Consequently the patterns 
and the associated localised eigenfunctions occur with a non-zero density along the 
infinite cluster and thus produce jumps of the IDS at the corresponding energy. In 
52 it is shown by physical arguments that the discontinuities of N constitute a 
dense set of energies. 

Actually, uniqueness of the infinite cluster is not used in the arguments of |52| 
and a similar argument for constructing finitely supported eigenfunctions does work 
on the Bethe lattice as well, although there the infinite cluster is not unique. For the 
quantum percolation model on the Bethe lattice locally supported eigenfunctions 
have been constructed in 06 § 7]. 

A mathematically rigorous study of the quantum percolation model on amenable 
graphs is undertaken in |314| . There the discontinuities of the IDS are explained 
in terms of the breakdown of the unique continuation property of eigenfunctions of 
the adjacency operator, see also Remark |3. 1.21 Moreover, the set of these energies 
is characterised in the case X = Z d . From a wider perspective, the properties of 
this set are related to the Atiyah conjecture, cf. |78) . 

Remark 3.1.4. While the continuity of the IDS has clearly to do with the 
distribution of eigenvalues of the random Hamiltonian, it only captures a part of 
the properties of this distribution. The theory of level statistics is concerned with 
the finer structure of the fluctuations of eigenvalues. It can be studied by an ap- 
propriate scaling procedure. This has been carried out for certain one-dimensional 
and discrete models in |228l 12271 126T1 12601 12251 H23] . 

3.2. Application to Anderson localisation. In the last paragraph the im- 
plications of Wegner estimates for the IDS were presented. Now we focus on the 
second main application of those bounds, namely Anderson localisation. 

As we discussed earlier in § 11.31 this phenomenon tells us that a random family 
of Schrodinger operators exhibits in a certain energy interval dense pure point 
spectrum, almost surely. Moreover, the eigenfunctions of the eigenvalues lying 
in this interval decay exponentially. Even a stronger property, namely dynamical 
localisation, holds. See § 11.31 for more details and references. 
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For multi-dimensional configuration space there are two techniques at disposal 
to prove localisation: the fractional moment method and the multiscale analysis 
(MSA). The first one is also called Aizenman-Molchanov technique and was intro- 
duced in |9"1 14[ |5 [ 11291 17]. It was so far applicable only to lattice Hamiltonians, up 
to the recent extension to continuum configuration space [6]. For discrete models it 
has in fact been proven |10L 111) that in the energy regime where the MSA applies, 
the Aizenman-Molchanov method works, too. 

However, we will discuss in a little more detail only the MSA, since it has found 
applications to a variety of models and since the Wegner estimate is a key ingredient 
in the MSA. We first sketch the basic ideas of the MSA, and then discuss shortly 
its history. 

To carry through the MSA one needs two a priori estimates: the initial scale 
estimate and the Wegner estimate. These two conditions essentially determine for 
which single site potentials u, single site distribution measures \i and which energy 
intervals localisation can be derived. Note that u and /i are parameters which 
determine our alloy type potential, see Definition If .2.11 

In the literature one can find multiscale analyses which are adapted to operators 
describing the propagation of classical waves or to abstract families of differential 
operators, see among others [TUH ITU51 I5T1 1571 12351 [T3UI IT77I . In this con- 
text one has also to make sure that certain other conditions are satisfied, like the 
geometric resolvent inequality, the generalised eigenfunction expansion, a rough up- 
per estimate on the trace of spectral projections of finite box operators (obtained 
e.g. by the Weyl asymptotics), etc. However, since we discuss here only (random) 
Schrodinger operators, these conditions are automatically satisfied, cf. (J95, § 3.2] 
or [TT9l App. A]. 

The multiscale analysis is an induction argument over a sequence of increasing 
length scales Ik, k £ N. Each scale lk+i is a power 1% of the preceding one, where 
a e]l,2[. Actually, for technical reasons one truncates the scales so that all Ik lie 
in 6N. 

One considers the restriction of the alloy type model H u to the open cube 
\( k ) ~ A; fc (0) of side length Ik- The corresponding restricted operator is denoted 

(k) 

by Hi, , where Dirichlct, Neumann or periodic boundary conditions ensure its 
sclfadjointness. One wants to study decay properties of the Green's function of 
Hu°\ i.e. the integral kernel of the resolvent operator R^\z) — (Hu ~ z)^ 1 , 
where z is taken from the resolvent set C \ a(H w k '). Since we are not interested in 
pointwise properties of the kernel of R u <(z), and since they tend to be unpleasant 
near the diagonal, we may consider instead the sandwiched resolvent 

X^iHl-E)- 1 ^ 

Here x° ut denotes the characteristic function of the boundary belt A;_i \ A;_3, and 
X m the characteristic function of the interior box A ; / 3 . 

The initial scale estimate is stated in terms of the notion of regular cubes. A 
box A; = A/(0) is called (E,j) -regular if I e 6N, E £ cr(^), and 

(3.6) ||x out (^-^)- 1 X in ||<e-' 1 '' 

The exponent has to satisfy 7 > W for some (3 > — 1. So regularity describes 
quantitatively how fast the Green's function on a finite box decays. The exponent 
(3 must be greater than —1 otherwise the rhs of l|3.6(l could be just one. 
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The initial scale estimate is satisfied if there exist a scale h < oo such that for 
some £ > 

(3.7) PH VE e I : A t is (E, 7 )-regular for oj} > 1 - Z -s 
for any l>l\. 

The weofc /orm of the Wegner estimate as it is needed for the MSA is: 

(3.8) ¥{d(a{H l J,E) < e~ l °} < Z~ 9 for all Z G 6N 

where 6 < 1/2 and q > d. 

The initials scale estimate (|3.7|l serves as the induction anchor of the MSA. 
The induction step uses the Wegner estimate and proves that the exponential de- 
cay property holds on the subsequent scale I2 with even higher probability, and 
that the decay exponent 7 essentially does not change. As one repeats the pro- 
cedure on the scales h,ls,... one proves that the decay of the Green's functions 
X° ut (Hu — E)~ 1 \ m holds with probability which tends to one, with error bounded 
polynomially in 17 1 . 

Thus one establishes the exponential decay of the sandwiched resolvent on arbi- 
trary large cubes, where the decay exponent 7 is bounded away from zero uniformly 
in the scales. Now one uses polynomial bounds on the growth of eigenfunctions and 
subsolution estimates to prove spectral localisation, cf. for instance 295, § 3.3]. To 
prove dynamical localisation one has to do more work, see e.g. [2951 § 3.4] or 

[TTsl 16711720]. 

The assumptions for the MSA depend on several parameters, and so do the 
various versions of localisation which may be obtained by it. Recently Germinet 
and Klein showed in |120| how to optimize the dependence of the MSA on the 
various parameters, i.e. how to obtain with the weakest assumptions in the input 
the strongest conclusions. 

The MSA was introduced by Frohlich and Spencer in [113] . The method 
applied to the Anderson model on the lattice and experienced various improvements 
and applications [2181 12191 ITT21 12831 173] , 

Based on results from [319] and [291] von Dreifus and Klein presented in [320] 
a streamlined version of the MSA. Although results on localisation for continuum 
Hamiltonians existed earlier [21711193] . it was this simplification of the MSA, which 
made alloy type Schrodinger operators more accessible to systematic research. 

There was a series of articles which proved various variants of the MSA for 
continuum models [TTHl l54l Il80l fTTol [TTHl Il09l [671 12951 Il20] . Other works 
concentrated onto identifying energy/disorder regimes where one can prove the 
assumptions needed for the MS A to work [T78l 11791 Il80l l23l ITTol IT691 Il20l 
294, 15H)1 15171 IT3H 15531 1153] . 

Remark 3.2.1. One dimensional models play a special role in the game of 
localisation. Namely, for d = 1 there are some specific techniques available which 
do not exist in higher dimensions, or are not as effective. Some examples are: the 
transfer-matrix method, study of the Ljapunov exponent, ODE techniques, Priifer 
transformation. 

Consequently, in one space dimension localisation has been proven for some 
models even before the MSA technique was available. See [1281 12281 [1271 liil 12631 

193 for various models on Z or M. Furthermore there are certain models which even 
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now can be treated only in one dimension. This applies to the following random 
Schrodinger operators: random displacement model |43LI284j . potentials generated 
by a Poissonian field |43| . alloy type potentials with changing sign (at all energies) 
|296j , discrete and continuous models with Bernoulli disorder |46l 12751 166) . This 
restriction to d = 1 is partially due to the fact, that there is no Wegner estimate at 
disposal in these cases. 

The following paragraph gives an idea where the Wegner estimate is used in 
the MSA. 

3.3. Resonances of Hamiltonians on disjoint regions. Rather than de- 
scribing precisely how the Wegner estimate enters in the induction step of the MSA 
we will confine ourselves to present an illuminative example due to Spencer 290 . 
It was originally formulated for lattice Hamiltonians, but can also be considered in 
the continuum case, as we have learned from P. Mvillcr. 

As we mentioned earlier (|3.2() implies for < S < 1 

(3.9) PH d(a(Hl), E)<6}< C W {E) S a |A,| fc 

This inequality implies that with respect to the parameter u> the eigenvalues of 
do not cluster on the energy axis. To give a more precise meaning to this 
statement, consider two box Hamiltonians H 1 = H^'^ and H 2 = \ Assume 
that the boxes A; (a;) and Ai(x) are sufficiently far apart such that H 1 and H 2 are 
independent. Let / be a bounded interval and consider the event 

where Cj stands for cr(H l ) n /. Let Ai, . . . , A^v be the eigenvalues of H 1 in /. By 
Weyl's law we know that N < Cj|Aj|, where Cj is independent of oj. Since 

N 

f2(<7i,cr 2 ) C [J £l(n,a 2 ) where Q,(n,a 2 ) := {u>\ B 2S {\ n (uj)) n a 2 ^ 0} 

n=l 

we may use Q3.9JI to conclude 

(3.10) P{n(<7i,«7i)} < CT Q |A / | 6+1 

This means that resonances of H 1 and H 2 , i.e. the occurrence of approximately the 
same eigenvalues for both operators are very unlikely. 

The feature which is common to Spencers example and the MSA is the effect 
of resonances between Hamiltonians which are localised to disjoint cubes. As we 
mentioned earlier, in the induction step of the MSA one puts together boxes A/ of 
side length I to form a larger cube A^ of side length L. Assume that one knows 
already that the Green's functions of the operators living on any one of the 
small cubes A; decays exponentially. 

The Schrodinger operator H% on A^ is obtained when we remove the boundary 
conditions which separate the smaller boxes A;. The question is whether the Green's 
function on Al will still decay (with approximately the same rate). To answer 
this question affirmatively it is not enough to know the exponential decay of the 
individual Green's functions on the small boxes A;, but it has to be ensured that 
they are not in resonance with each other. 
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Resonance means in this context that the spectra of two restriction H 1 , H 2 of 
Hu to disjoint cubes are very close to each other, and can be formulated quantita- 
tively in terms of 

(3.11) d(a 1 ,<7 2 ) := inf{d(Ai, A 2 )| Ai G (T 1 ,X 2 E cr 2 } 

as we will see below. 

The model situation we are about to consider is easier than the one occurring 
in the MSA because we do not introduce boundary conditions but confine ourselves 
to the analysis of the discrete spectrum below zero. 

Example 3.3.1 (Spencer's example |290l p. 903-904]). Consider two smooth 
potentials V\ , V 2 < with compact support and set 

(3.12) Vi := suppV* C B r (ai),r > 0,a t e R d ,i = 1,2 

It follows d(Vi, V2) > \a\ — a 2 \ — 2r =: g. Consider furthermore the operators 

(3.13) H := H + V 1 + V 2 , H := -A 

(3.14) Hi:=-A + Vi, i = 1, 2 

Denote by Oi := cr(Hi) fl] — 00, 0[, i = 1,2 the negative spectra, which are 
purely discrete. We are interested in the localisation and decay properties of the 
corresponding eigenfunctions. 

We look at two cases, where the first has a special symmetry and the second 
corresponds to the situation one expects to occur in a random medium. 

Case (A): 

Consider first the exceptional case in which V 2 is obtained from V± by a reflection 
along an axis of symmetry. Without loss of generality 

(3.15) V 2 (x 1 ,x 2 , . . . ,x d ) = Vi{-x\,x-2,...,x d ) 
Thus H commutes with the reflection operator 

(3.16) n: L 2 {R 2 ) ^L 2 (M 2 ), {Uf){x u x 2 , . . . , x d ) = f{-x u x 2 , . . . , x d ) 
In particular, for every eigenfunction -0 

Hip = Xip, A < 

the reflected function Hip is an eigenvector of H as well. If ip is localised around ai 
Hip will be localised around a 2 . Thus a typical vector form spanj^, Hip} will have 
non negligible amplitudes both at a\ and a 2 , even for large distances g. For short, 
eigenfunctions of H do not need to have just one centre of localisation. 

We are dealing with a resonance between the two disjoint regions Vi,V 2 , or 
more precisely between the spectra of Hi and H 2 . Actually, we encountered the 
extreme case where the the spectra <j\ and a 2 are not only close to each other but 
identical. 

The example we just considered exhibited a special symmetry, namely [H, II] = 0. 
For random potentials we expect generically that such symmetries are absent and 
that the spectra o\ and o~ 2 have positive distance. This situation is considered in 
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Case (B): 

We give a condition on d(o"i,<72) which ensures that the eigenfunctions of H (de- 
fined in (|3.13f) ') are localised at only one of the potential wells Vi,V2- Namely, 
assume that 

(3.17) d{a 1 ,a 2 )>e-^=:e 

For an eigenvalue A < — g^ 1 , with corresponding equation Hip = Xip, we have either 

(3.18) \X{ - A| > e/2, VA{ € CTl or |A 2 ' - A| > e/2, VA J 2 € ct 2 
Assume without loss of generality the first case. The eigenfunction equation implies 

-ill = (H, - A)- 1 ^ 
Applying twice the resolvent equation we obtain 

(3.19) -i\> = [i?o - RoViRo + RoV.R^Ro] V 2 ip, 

where Ri := [Hi — A) -1 , i = 0, 1 denotes the resolvents. We show that the am- 
plitude of ip on Vi is exponentially small in the parameter g. Denote with x* the 
characteristic function of Vi for i = 1, 2 and multiply (|3.19|l with x 1 

(3.20) - xV = X^oX 2 !^ - X^o^ix^oX 2 ^^ + x^oV^V^RoX^ip. 
The free resolvent decays exponentially, see e.g. [3] or [2551 IX. 30], 

R 1 ' 2 :=\tfRoX , \\<e-°S=X, 
the terms 11^2^11 , \\x 1 RoVi\\ are bounded uniformly in g, and (|3.18(l implies 

H-Rill < \ =2e^ 

Consequently 

HxVll < R 1 ' 2 \\V 2 *l>\\ + W^-RoVtW R 1 ' 2 \\V 2 ip\\ + Hx^o^ill HfliViH iZ 1 ' 2 \\V 2 ip\\ 
is bounded by a constant times exp(— ^fg) \\ip\\, since g^ 1 < —A. 

Let us finish this section by discussing some aspects and contrasts of the two 
cases considered in the example. 

(i) In general, the spectrum alone describes only general properties of the solution 
of the eigenvalue equation. In our example it is the additional information 
contained in H3.15|l and H3.17|l . respectively, which allows us to analyse the 
eigenfunctions more precisely. 

(ii) Obviously, in Case (A), the Green's function decays in space, too. However, 
this decay is not yet felt at the scale g, since \ip(ai)i(}(a 2 )\ converges to a 
positive constant for g — ► oo. On the contrary, in Case (B), the amplitude is 
small either at a\ or a 2 , so 

(iii) A semiclassical analysis of double well potentials is carried out, for instance, 
in [135] . 
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4. Wegner's original idea. Rigorous implementation 

In this section we present a proof of Wegner's estimate following his original 
ideas in [325 . His proof was originally formulated for the discrete Anderson model. 
In the meantime, it has been cast into mathematically rigorous form and adapted 
for continuum Hamiltonians. We follow mostly the arguments of Kirsch |159j . 
There are proofs of Wegner's estimate by other authors, which make use of the 
ideas in |325] . Let us mention [2Tfl 12161 HTSl ITT91 Il08l 1571156]. 

The theorem to be proven is 

Theorem 4.0.1. Let be as in Definition \1.2.1\ and assume additionally that 
there exists an k > such that 

(4-1) U > KX[_l/2,l/2]<* 

Then for all Eq € R there exists a constant Cw = Cw{Eq) such that for all I e N, 
E < E Q and all e £ [0, 1] 

(4.2) E{Tr[PU{E-e,E + e})]} < C w e l 2d 

This theorem is proven in the next section. Its bound with respect to the 
volume term l d is quadratic and does not yield a continuity statement for the IDS. 
Subsequently we show how this estimate was improved in |6U| . Denote by lo + and 
c<j_ the largest, respectively the smallest value a coupling constant may take. 

4.1. Spectral averaging of the trace of the spectral projection. We 

show that the expectation over the randomness smears out the eigenvalues of H l u 
and thus regularises the trace P^ (I) . 

By definition the spectral projection P^(I) = Xi{H l w ) is the characteristic func- 
tion of Hlj . For certain purposes it will be necessary to differentiate this function 
with respect to the energy parameter, which motivates the introduction of the fol- 
lowing smooth 'switch function'. 

Let p be a smooth, non-decreasing function such that on ] — oo, — e] it is iden- 
tically equal to —1, on [e, oo[ it is identically equal to zero and. 1 1 // 1 1 oo 

< l/e. Then 

XiE- e ,E+e[(x) <p(x-E + 2e) -p(x-E-2e)= dtp'(x-E + t) 

J-2e 

Thus by the spectral theorem 

PU]E-e,E + e[)< dt p' (H l u - E + t) 

J-2a 



in the sense of quadratic forms. Since B e (E) =]E — e, E + e[ is bounded and o-(H^) 
discrete, the above operators are trace class and we may estimate: 



Tr Pl(B t {E))] < Trf / dt p' (H l u - E + 1)] = V / dtp'(\ l n (u) -E + t) 

J -2t J „^R,T»'-2e 



" 2e n&l J - 2t 



where Ajj(a;) denotes the eigenvalues of enumerated in non-decreasing order 
and counting multiplicities. Only a finite number of terms in the sum are non-zero. 
More generally the above arguments prove the following 
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Lemma 4.1.1. Let H be an operator with purely discrete spectrum. Denote by 
Ai < A2 < • • • the eigenvalues of H. Then for £gK and e > 

Tr 



XB t (E){H)] <E / dtp'(X n -E + t) 



In the following we analyse the behaviour of the spectrum of the Schrodinger 
operator under the perturbation ujj u(- —j). Fix a box-size I € N, a lattice site j e A 
and a configuration of coupling constants w e O and consider the one-parameter 
family of operators 

(h H t := H + tU, where H = H l u and U = u(- - j) 

By the arguments in Ml. 21 the single site potential is infmitesimally bounded with 
respect to H, thus H t forms a holomorphic family of type (A) in the sense of Kato 
|152 | for t in a neighbourhood of the real line, cf. e.g. XII. § 2 in |256| . Moreover, 
H t has compact resolvent by XII. § 14 in 25£>]. Hence one may apply a theorem of 
Rellich 259 , see also Theorem VII. §3. 9 in 1 52 . It says that the eigenvalues and 
eigenvectors of H t can be chosen to be real analytic on R. Actually, each eigenvalue 
is holomorphic on a neighbourhood of R in the complex plane, but their intersection 
may contain only R. 

If X n (t) is a non-degenerate eigenvalue of H t , first order perturbation theory 
tells us that there exists a normalised eigenfunction ip n (t) such that 

(4.3) ^(to) = MM, Uip n (t )} 

Remark 4.1.2. This is sometimes called Hellmann-Feynman formula, and it 
holds true also if the eigenvalue A„ happens to be degenerate at t = to, cf. for 
instance |146| . One has however to chose the enumeration of the eigenvalues X n 
and eigenvectors ip n in such a way that the pair X n (t) , ip n (t) , t < to continues 
holomorphically into X n (t), ip n (t), t > to- Note that this is actually not the case 
with the enumeration we chose earlier, where X n (t) denotes the n-th eigenvalue of 
H t . There are two possibilities to solve the problem: either one chooses a somewhat 
unintuitive enumeration of eigenvalues which makes them — together with the 
eigenvectors — holomorphic functions of t. Or one sums over the eigenvalues. 
Namely, formula H4.3(l remains true if we sum over all eigenvalues which correspond 
to a degeneracy. More precisely, for a degenerate eigenvalue X n (to) denote by 
l,i;eN the largest numbers such that X n -i(to) = ■ ■ ■ = X n (to) = ■ ■ ■ = X n+ k(to) 
and set S(t) = Y^m= n -l ^m(*)< Then we have 

dS nJrk 

-(to) = W'mCfo)) Ulp m (t )) 

m—ri—l 



dt 



In the application of l|4.3|) in the next proposition we will be considering all eigen- 
values below a certain energy. Thus if we consider one eigenvalue participating in 
a degeneracy we will actually take into account all participating eigenvalues. 

The results for one parameter families of operators carry over to the multi- 
parameter family to 1— > . Thus we have 



^ du) j ^ 
jeA jeA 



(lpn,u(- -j)lp n ) 
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where ip n are normalised eigenvectors corresponding to A n (£/,",). By assumption 
j|4.1(l we have 

(4.4) X^r,, «(■-•?>«) > K >° 

j£A 

Now the chain rule 

^ = ^ , ( A »™-- B + t )2^-^- 

implies 

dp(X n (Hl)-E + t) 



(4.5) p'^HD-E + t) <^J2 



Due to monotonicity, integrating over one coupling constant we obtain 

dp(\ n (Hp-E + t) f dp(\ n (Hp-E + t) 

duj -WJWoo J u 3 

= ll/l|oo[p(A„(w, j = max) — E + t) - p(\ n (u,j = min) - E + t)] 

where \ n {ui,j = max) denotes the n-th eigenvalue of the operator 

H l u (j = max) := H l u + (lu+ - Wj ) u(x - j) 

where ujj takes its maximal value. Analogously we use the notation \ n {ui, j = min) 
This proves 

Proposition 4.1.3. 
E{Tr[P<([£-e,£])]} 
"/II. 



< 



— Y\ I dt J2 E{p[X n (u),j = mzx)-E+t]-p[\ n (iJ,j = win)-E+t]} 



The upper bound can be also written as 

l/l' ' 2 



^ [ dt V E f Tr [p[Hl (j = max) -E + t]-p[Hl (j = min) -E + t]]} 

K J-2e t ' 

Since p < 



^ p[A„(w, j = max) -£? + *]- p[A„(o;, j = min) — E + t] 

riGN 

(4.6) < - ^ p[An(w, j = min) -£? + *]< C B+3 e Z d < C' Eo+3 l d 

by bound (hi) in Lemma \l. 7. 31 This proves Theorem l4.U.ll 

Remark 4.1.4. The suboptimality of the volume bound in Theorem 14.0.11 is 
due to the rough estimate (|4.6|l . The right hand side of the inequality is the net 
increase of the number of eigenvalues in the energy interval ]E — t — e, E — t + e[ due 
to the decrease of the j-th coupling constant from its maximal to its minimal value. 
This quantity is expected to be independent of A. However, in 1)4. 6J) we estimated 
it by the total number of eigenvalues below the energy E + 3e, which is by Weyl's 
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law proportional to the volume of A. Thus we get an extra volume factor in the 
upper bound of the Wegner estimate. 

4.2. Improved volume estimate. In 160] Combes, Hislop and Nakamura 
obtained a Wegner estimate analogous to the one in the last section, but with 
upper bound linear in the volume. The Wegner estimate they proof is somewhat 
weaker in the energy parameter. 

More precisely the main result in 60 may be formulated as follows. Let 1q := 
]E-,E+[,Ig H <j(Hq) = be a spectral gap of Hq. Denote by H^ the operator 
Ho + X^fceA, LJ fcw(- — k) and by the corresponding spectral projector. 

Theorem 4.2.1. Let H u be as in Defjnition \1.2.1\ and assume additionally that 
the periodic potential V pcl - is bounded below, < u £ C C (1R ) and u is not identically 
equal to zero. Let E £ Ig, < eo •= \d{E, Iq) and a < 1 then there exists a finite 
constant C such that 

(4.7) E{Tr[Pl(B e (E))]} < C e a l d for all < e < e and I £ N 

Note that H l u is not an operator restricted to a bounded domain. This is 
the reason, why the theorem does not apply to energies inside the unperturbed 
spectrum a (Ho). Estimate (|4.7|l implies the Holder continuity of the IDS outside 
the spectrum of the unperturbed operator Ho- A sharper bound, linear in e, was 
obtained in [23] . However, its proof is based on a different technique, cf. Section^ 
We will prove here a slightly simpler fact than the above, namely 

Theorem 4.2.2. Let H^ be as in Definition \1.2.1\ and assume additionally that 
the periodic potential V pcl is bounded below, u £ L^(M rf ),p = p(d) and for some 
k > 

w > K X[o,i] d 

Then, there exists for any Eo £ M and a < 1 a finite constant C such that 

E {Tr[P^(B e (E)]} <Ce a l d 
for all e £ [0, 1], E < E and I £ N. 

This result is proven in Section 3 of 189 using partially different methods. 
Note that the earlier papers |193l 151] contain sharper estimates, cf. Section[S] 

The fundamental contribution of |60| was that it replaced the Weyl-type volume 
estimate (|4.6ll . whose drawback we explained in Remark 14.1.41 We present now the 
improved volume estimate of 60.. Together with the argument from the previous 
§ 14. II it proves Theorem [4 .2. 21 

The trace 

(4.8) Tr[p[H l u (w,j = max) -E + t]- p[Hl(w,j - min) -E + t}] 

can be expressed using the spectral shift function £ = £(•, H + U, H), abbreviated 
SSF, of the pair of operators 

H := H u (oj,j — min) and H + U, where U — — oj_) u(x — j) 

The necessary estimates on the SSF are collected in Appendix [S] Since the differ- 
ence of our operators is not trace class, we have to use the indirect definition of the 
SSF by the invariance principle, cf. I]A.9]) . Let Co £ M be such that iJ w ,i7o > Co 
for all uj and g(x) = (x — Co + l)~ k for some k £ N. 
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For k > ^^,k e N the operator g(H + U) — g{H) is trace class and the 
invariance principle implies 

Tr(p(H + U — E — t) — p(H — E — t)) = — f p'(X)((g(X),g(H + U),g(H)) dX 

(4.9) <-(jp'{XydX) 1,q (f t(g(X),g(H + U),g(H)) P dx) 1/P 

In the last line we used the Holder inequality and p,q > 1 are conjugate exponents 
i + i = 1. Remember that we choose p depending on e. Thus, its derivative is 

bounded by e _1 times a constant and 

(fp'(XydX) 1/q < C(l/e)^(Jp'(X)dX) 1/q = Ceh 1 = C ^ 

Since on the support of p' the function g is uniformly bounded away from zero a 
transformation of variables gives 



U//< A).//. // + r ).,/.//) >" dX)'" <c( [ C(A, g(H + U), g{H)f dX ' ' 



suppp 7 



= C \\g(H + U)-g(H)\\ 1 / i P /p <C 



where we used in the last line Theorems IA.1I and IA.6I from the Appendix. There 
one can find the definition of the super-trace class ideal J\/ p and its norm. The 
constant C is independent of the box A/, the lattice site j, the configuration of the 
coupling constants w^, k ^ j, and of e. 

Hence we have a volume independent bound for (|4.8|) . in contrast to the es- 
timate discussed in Remark 14.1.41 The bound is Holder continuous in the energy 
parameter. In view of Proposition 14. 1 . 3l the proof of Theorem l4.2.1l is finished. 



Remark 4.2.3. Recently Hundertmark, Killip and the author |137j found a 
different, shorter way to prove the super-trace class estimates and to apply them to 
bound the SSF. The basic observation is, that one can control the singular values 
of the difference of two Schrodinger semigroups. In fact the singular values decay 
almost exponentially, and the semigroup difference is therefore in any super-trace 
class ideal. 

We state without proofs their result on the decay of the singular values and the 
estimate on the SSF it implies. For simplicity we consider here the case where the 
magnetic vector potential is absent. 

Theorem 4.2.4. Let Hi = -A+V andH 2 = Hi+u, withV,u e Ll oc (R d ),V,u > 
— ^Cq. Denote by H\,H\ the corresponding Dirichlet restrictions to the cube A/ = 
[—1/2, l/2] d . Set V^a :— e~ H ± — e~ H i . There are finite positive constants C\, c 2 such 
that the singular values p n of the operator V^ ff obey 

(4.10) Pn<C ie - C2nl/d 

The constants depend only on d, Co and the diameter of the support of u. 
Let p be a switch function as above. 
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Theorem 4.2.5. There is a constant depending only on d, Co, diamsuppu and 
E + e such that 

Tr [p(H l 2 -E)- p(H{ - E)] < const log(l + l/e) d 

for all e > and I £ N. 

4.3. Sparse potentials. Form the physical point of view there are some 
interesting models which have a potential 

(4.11) V u (x) = J2 u ku(x-k) 

ker 

resembling the alloy type model. However, the set T may be much more general than 
the lattice Z d . A class of particular interest are surface models where T = {0} x % v 
and v < d is the dimension of a hyperplane in whose neighbourhood the potential 
is concentrated. The literature on such models includes |48L 11881 11891 136] . see 
also [961 19 7| . The results in this paragraph are taken from |172| . 

Here we will consider arbitrary sets T, which are uniformly discrete in the 
following sense 

sup #{r n B x (x)} < oo 

xeR d 

For uniformly discrete T the number of points of T contained in the cube A; (a;) can 
be bounded linearly in the volume of the cube and independently of its centre x. 

Consider a background Schrodinger operator H = — A + Vp Cr with a periodic 
potential Vp Cr € L^ nit loc (K d ) where p = p(d) is as in (|1.1|) . By adding a constant 
we may assume that inf o~(Hq) = 0. Let H u = Hq + V u be an random operator 
with an generalised alloy type potential as in (|4.11|) . 

As before stands for the restriction of H u to the cube A; with Dirich- 
let boundary conditions (we may as well use Neumann or periodic ones), and 
denotes the corresponding spectral projection. 

Theorem 4.3.1. Assume that the single site potential u £ L^(M. d ) is non 
positive, and that the single site distribution /i has a density f £ L£°([0, oo[). 

Then, for any a < 1 and ~ E' < there exists a finite C such that for any 
E £ M, e > satisfying E + 3e < -E' : 

E [TrP^(B e (E))] <Ce a l d 

The proof of the theorem follows from the arguments of § 14.11 and § 14.21 once 
a replacement for the estimate (|4.5f) has been established. This is provided by the 
following lemma. Denote by A^ = {k £ T\ suppu(- — x) D A; ^ 0} the set of indices 
whose coupling constants influence the value of the potential in the cube A; . Recall 
that the supremum of the support of / is denoted by lu + . 

Lemma 4.3.2. Assume that the n-th eigenvalue of the operator H l u satisfies 
X l n (oj) < -E' < 0. Then 

X- dp(\ l n (u;)-E + t)' 



pW^-E + t)^ 1 ^ 



, AO. ^IVk 

feeA+ 



Proof. Let ip n be the normalised eigenfunction corresponding to X l n (u>). Then 
V>n satisfies by definition {ijj ni (Hq — \ l n (uj))ip n ) = — {i/; n , V w ip n ). We have 

^ ^fc (i>n,-Uk(- ~ k)lp n ) = -(ipn, V^lpn) = (lp n , (H l Q - \ l n (uj))ip n ) > E' 
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Now we have by the Hcllmann-Feynman theorem 

dx l n (cj) 



E' 



V" = ^ (^n,-Uk(—k)lp n ) > V" UJk (lpn,-Uk( — k)lp n ) > 

* — ' OUJl. * — ' A — ' w + 
fcSA+ fceA+ keh+ T 



This gives 

p'(\ l n (u)-E + t) = 

(4.12) < 



E 



L fceA+ 



duj k 



E 

fceA+ 



dp(\ l n {uo)-E + t) 



duj k 



E' 



E 

fceA+ 



dp(\ l n (v)-E + t) 
duj k 



□ 



Note that since p is monotone increasing and u is non-positive, (|4.12() is a 
non-negative real. 

Related models and results as presented in this paragraph are discussed in 
Section 3.1 of |57| . In |172j two more classes of generalised alloy type models are 
analysed. Firstly, the case where the number of points in T n B\{x) is not uniformly 
bounded, but grows at a controlled rate as x goes to infinity. Secondly, the case 
where Y is itself a random point process, for example of Poissonian type, cf. also 

m 

4.4. Locally continuous coupling constants. In this paragraph we present 
a Wegner estimate which requires the coupling constants ujk to have a continuous 
distribution merely in a neighbourhood of their extremal value uj + = sup supp /. 
Both the result and its proof are taken from |172| . 

Consider a background Schrodinger operator Hq = — A + V pcr with a periodic 
potential V por £ LjJ nif loc (]R <i ). Let H u = H + V u be an random operator with an 
alloy type potential V w . Assume that the coupling constants Uk, k £ Z d take values 
in the bounded interval [w_,cl> + ]. By modifying the periodic background potential 
we may consider only the case that the coupling constants are non-negative. For a 
value lu c £ [0, w+] introduce the auxiliary periodic potential V c = uj c J2kez. d u ( x ~k) 
and the threshold energy E c — inf a(H + V c ). 

Theorem 4.4.1. Assume that the single site potential u £ L£°(R ) is non 
positive, and that the restriction of the single site distribution p c := has a 

density f £ L°° . 

Then, for any a < 1 and E' < E c there exists a C such that for any E £ R, e > 
satisfying E + 3e < E' : 

E [Tr P^(B e (E))] <Ce a l d 

Proof. The value u> c is a critical one for the random variable ujk in the sense 
that for ujk > uj c we know that it is continuously distributed, while for smaller 
values we do not know anything. We introduce a corresponding decomposition of 
the 'probability' space := x feey \+R — R L - This is the part of the randomness on 
which the restricted Hamiltonian depends. For a given configuration of coupling 
constants {^k}kGA+ se t 

A ac (uj) ={k£A+\uj k > lu c ] 
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This defines an equivalence relation on (li by setting for any A C A + 

(1(A) := H A ac (uj) = A} 

Consequently 

(4.13) Yl f L II d Vk(uk)Xn(A)(u>) = 1. 

Split the potential now into two parts, a singular and an absolutely continuous one. 
The singular one 

v u( x ) : = X! w fc ufc(a;-fc)+ 2J u c u k {x - k) > V c (x) 

will be considered as part of the background operator, while the absolutely contin- 
uous one 

V£°(x) ■= r k u k (x-k) = ^ r k u k (x~k), with r k = uj k - qj c > 

will be used for spectral averaging. 

Consider an eigenvalue X l n < E' < E c and an eigenfunction H^ip n = X l n ip n and 
set 6 = E C -E'. We have 

-<</>„, KT</>„) = (#o + V£ - A^„) > (^„, (H l + V C - \ l M > ,5 

which implies similarly as in Lemma 14.3.21 

Consider first the case A + and estimate 



/ LI ^(^fckn^M^ / p'(K(w) - E + t) 

jRL kEA+ nGN 

< ^? — - / FT rf i"h)Xa(A)M)] / 



dp(X n (uj) - E + t) 

jgA ac J 



As we know that all sites j S A ac correspond to coupling constants u>j with values 
in the absolutely continuous region of the conditional density / we may estimate 
as in § O 

- 2^ / d MK)xn(A)(w) — 
...Ji 



vie! 



, 0p(X„(w)-.E + i) 



= -£/ + /(^-)^ 

< || /|| oo £ [p(A„(w,tJj = u> c ) - E + t) - p({\ n (uj,ujj = u> + ) - E + t)] 



ufcf 



which can be estimated as in § 14.21 We have to say something how we deal with the 
special case A = 0. In this situation V^ c = and H u = H + V* > H + V c > E c . 
Thus there are no eigenvalues in the considered energy interval for this potential 
configuration. 
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Finally we use the decomposition l|4.13|) to finish the proof: 
E (TrPU[E -e,E + e])) 

< E / II M"k)Xn(AMj2 iW+ ~ Wc) te\\f\\ooC(a)e 



a-l 



AcA+"" keA+ jeA 
< 4G(a) llJllooE * 



□ 



4.5. Potentials with small support. In § 14.11 we used in a crucial step 
in the derivation of the Wegner estimate that the single site potentials were lower 
bounded by a partition of unity 

^2 u(x-k)> k on R d 

feez d 

It is of natural interest, whether a Wegner estimate holds if this condition is relaxed. 
In this paragraph we consider the case that u is of fixed sign, but has small support. 
More precisely, we assume throughout this paragraph merely that there is an open 
set O C K d and a positive k such that 

(4.14) u(x) > k X o 

The first Wegner estimates under this relaxed condition on the single site potential 
were derived for spectral boundaries, i.e. for energies either near the bottom of the 
spectrum, or near an internal spectral boundary. The case of the infimum of the 
spectrum was treated e.g. in [18011159] . and internal spectral boundaries in |170j . 
These works derived a Wegner estimate where the volume dependence of the bound 
was growing faster than linearly. Thus they were not sufficient to derive a result on 
the regularity of the IDS, cf. our discussion in § 13.11 A linear bound for the same 
energy regimes was found in |23l I60| . 

By now there are Wegner estimates which under the relaxed condition (|4.14(l 
derive Wegner estimates valid for any bounded interval on the energy axis. Wc 
consider first the one-dimensional case where the result is particularly clear and 
the proof simple. We follow |171| in the presentation, see 122, 56 for other 
proofs. 

Assume that the single site potential u and the periodic potential V^er are 
bounded. 

Theorem 4.5.1. Assume that u is compactly supported and obeys l|4.14|l . For 
any there exist a constant C such that 

(4.15) E [Ti P^{B e (E))] <C el, Ve G [0, 1},E < E ,l £ N 
Thus the IDS is Lipschitz-continuous. 

Proof. First we show how to replace (|4.4() in the case of small support. By 
shifting the origin of M. d we may assume without loss of generality that there is a 
s > such that A s (0) C O. Likewise, we may assume k = 1 by rescaling the single 
site potential and the coupling constants. 

We set S — UfcgA A s (fc). The Hellmann-Feynman theorem gives us 

keA " fceA 
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where tp n is a normalised eigenfunction corresponding to X l n (oj). 

If the integral on the right hand side would extend over the whole of Aj it would 
be equal to 1 due to the normalisation of tp n . A priori the integral over 5* could be 
arbitrary close to zero, but the following Lemma shows that this is not the case. 

Lemma 4.5.2. Let I be a bounded interval and s > 0. There exists a constant 
c > such that 



A s (fe) 



M 2 > c 



Ai(fc) 



for all I € N, all k € A; and for any eigenfunction ip corresponding to an eigenvalue 
Eel of Hi. 



Proof of the Lemma. For 



A 3 (k+y) 



dx\ip{x)\ 2 = 



A 3 (fe) 



dx \ip(x - y)\ A 



one has 
d 



dy 



dx 



A s (fc) 



d_ 

dy 



ip(x - y) 



ip(x -y) + 



8 



dxtp(x - y) -^-tp(x - y) 
A 3 (fc) cy 



< 



2 IMIi 2 (A s (fc+ a )) ll^'llL 2 (A s (fc+y)) ■ 

Sobolev norm estimates (e.g. Theorems 7.25 and 7.27 in [T24 ) imply 

||^'IU 3 (A s (fc+2/)) < C 5 W^\\mA s (k+y)) + ||V>"||i 2 (A s (fc+y)) 

By the eigenvalue equation we have 
d 



(4.16) 



dy 



4>{y) 



<c. 



6 H\\h { A s ( k+ y)) = Ce <Ky), C 6 = C 6 (||F per + V»- E\\ x ) 



Gronwall's Lemma implies (j>(y) < exp(C6|y|) (f)(0) and thus 



Ai(fc) 



|0| 2 <e 



c 6 8 -i 



A a (fe) 



a 



Thus J s | ■01 2 > c J A; 1^1 2 with the same constant as in Lemma T4.5. 21 
It remains to estimate the spectral shift 

(4.17) J2 [P^n^'i = max ) -E + t)- p(X l n (uj,j = min) -E + t)] 

We may assume without loss of generality that the single site potential u is sup- 
ported in the interval [—R,R]. Introduce now the operator Hl D (j = max) which 
coincides with (j — max) up to additional Dirichlet boundary conditions at the 
points j — R and j + R. Likewise, H^ N (j = min) coincides with H^j = min) up 
to additional Neumann boundary conditions at the same points. Their eigenvalues 
are A^ D (w,j = max) and X\l {ui,j — min), respectively. By Dirichlet-Neumann 
bracketing, the square brackets in i|4.17[) are bounded by 

(4.18) p()!f{u,j = max) -E + t)- P (X^ n (lu,j = min) -E + t) 
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Since for both * — D,N the Hamiltonian HJ' is a direct sum of an operator Hfc* 
acting on L 2 (j — R,j+R) and another one H%* acting on L 2 (A; \ \j—R,j+R\) the 
sum over the terms in (|4.18|) can be separated: 

(4.19) J2 P^n D ^) -E + t)- p{Kl N {u) -E + t) 

n 

(4.20) + Y,P( X n D (u>3=^)-E + t)-p{\^ N {u,j = vAn)-E + t) 

n 

Note that the eigenvalues in l|4.19[) are independent of uJju(- — j). Since the differ- 
ence in the boundary conditions is a rank two perturbation in resolvent sense (see 
e.g. |282| L the interlacing theorem says that 

p(A^ D (w, j = max) -E + t)< p(E c n ? 2 (u;,j = max) -E + t) 

A telescoping argument bounds the whole sum in l|4.19|l by twice the total variation 
of p, which is equal to one. The sum in 14.20fl we estimate by 

Tr [x[E-3e,oo[(Hj!j D {j = max)) - X]E+3e,oo[(H^ N (j = min))] 

< 2 + Tr [ X[ E-^oo[m D ^ = min ) + Moo) - X]E+3e,oo[( H i' D U = min ))] 

which is bounded by a constant, that is independent of A;, j £ A/ and e > (on 
which p depends). □ 

In the remainder of this section we give an overview of various Wegner estimates 
which are based or related to techniques presented at the beginning of Section 0] 
However, we refrain from giving the proofs of this results but refer to the original 
articles. 

In |56j Combes, Hislop and Klopp study multi-dimensional alloy type models 
with single site potentials of small support, and establish the Holder continuity of 
the IDS at all energies. They consider the case where the single site potential u £ 
L^°(E d ) is non-negative and not identically equal to zero, and treat three different 
situations. In all of them the unperturbed background operator Ho — (— zV — A) 2 + 
Vq may include a magnetic vector potential A and a (scalar) electric potential Vo- 
They have to satisfy some regularity conditions such that Hq is selfadjoint and has 
C^°(K. <i ) as an operator core. The coupling constants are distributed according to 
a bounded, compactly supported density. 

(i) The background operator Hq has an IDS No, which is Holder continuous 

\N (E 2 ) -N (£01 < C |£ 2 -£ir 

with Holder exponent a €]0, 1]. The constant Co = Cq(I) can be chosen 
uniformly for E2, E\ in a given compact interval /. 

(ii) The background operator is periodic with respect to the lattice Z d and has 
the unique continuation property, cf. for instance 327 . The set {x\ u(x) > 0} 
contains an open subset of M. d . 

(iii) Let the space dimension be d = 2. Let Ho = (— iV — A) 2 + V pcr be a Landau 
Hamiltonian with vector potential A(x\,X2) — ^{—X2,x\) where B > is 
the (constant) magnetic field strength. The magnetic flux trough a unit cell 
satisfies the rationality condition 

(4.21) Be2nQ 
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The scalar potential V pC v is a Z d -periodic function in Lf oc (R d ). 

In case (0) set a c = otherwise a c = 1. 

In a follow up work |58] on the Landau Hamiltonian in collaboration with 
Raikov condition l|4.21(l has been removed. 

Theorem 4.5.3. Let 77 w be an alloy type model satisfying either one of the 
above conditions lfi)l - lfm)l . Then, for each a e]0, a c [, the IDS of H u is Holder 
continuous at all energies, with Holder exponent a. 

4.6. Holder continuous coupling constants. There is special interest to 
extend the known Wegner estimates to coupling constants with singular distribu- 
tion. The reason is the intuitive interpretation of the coupling consonants as nuclear 
charge numbers modulating the strength of atomic potentials. In this case their 
distribution would correspond to a pure point measure. 

So far the best result in this direction for multi-dimensional alloy type models 
was obtained by Stollmann in 294 . In Remarks 13.1.11 and 13.2.11 we mentioned 
already results for one-dimensional models with singular randomness. 

Stollmann's result applies to a single site measure /i which has compact support 
and which is merely Holder continuous. For e > denote 

s(e) := s(p,e) := sup{/i([a, b])\ b - a < e} 

Theorem 4.6.1. Let H u be an alloy type model as in Definition \1.2.li but let 
the single site measure be merely Holder continuous. Assume additionally that the 
single site potential obeys u > X[o,i] d - Then for any E £ R there exists a constant 
C such that for any open interval 7 c] — 00, E[ and any I € 2N 

PHa(77i) n7^0}<C S (|7|)Z 2d 

holds. 

4.7. Single site potentials with changing sign. First Wegner estimates 
for indefinite alloy type potentials were derived in |180j . In |134| Hislop and Klopp 
combine the techniques from 180 and 60 to prove a Wegner estimate valid for 
general indefinite single site potentials and for energy intervals at edges of a(H u ). 
They assume the single site potential u 6 C c (IR rf ) satisfies u(0) 7^ 0. The density 
/ 6 of the random variable ujq (which may be in fact the conditional density 
with respect to lw 10 := (wfc)fe^o) is assumed to be piecewise absolutely continuous. 
For any a < 1 and any compact energy interval 7 strictly below the spectrum of 
the unperturbed operator 77o they prove 

P{a(77^)n7^0}<C|7| a / ti 

where the constant C depends only on a, d and the distance between the interval 
7 and er(77o). With a sufficiently small global coupling constant A the same result 
holds for the operator 77 + AVL for 7 in an internal spectral gap of 77 . The results 
of |134j extend to more general models including certain operators with random 
magnetic field. 

In § 15. 51 we discuss in more detail an alternative technique to obtain a Wegner 
estimate valid for single site potentials which change sign. It applies to a more 
restricted class of potentials but yields stronger results. In particular, it proves the 
Lipschitz continuity of the IDS at all energies. 
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4.8. Uniform Wegner estimates for long range potentials. Kirsch, Stoll- 
mann and Stolz proved in |169j a Wegner estimate for single site potentials which 
do not need to have compact support, but merely need to decay sufficiently fast. 
They consider u of polynomial decay 

(4.22) \u{x)\ < C(l + \x\ 2 )- m ' 2 

where m > is required to be larger than 3d. For certain applications they can 
also deal with the case where m is only larger than 2d, cf . [1681 1329) . 

For such single site potentials the restrictions of the alloy type potential to 
two finite cubes may be correlated, even if the cubes are far apart. This makes 
it necessary to use a enhanced version of the multiscale analysis for the proof of 
localisation. Among others, this requires a uniform Wegner estimate. By this we 
mean a Wegner estimate for the Hamiltonian H l u restricted to the cube A; which 
is uniform in the coupling constants ojk with index ||fc||oo > t where r is a function 
of I. 

To formulate the Wegner estimate from |169| let us first introduce some nota- 
tion. For any cube A C K d and A = A n 7L d we denote by IIa the projection 

n A : fi i-» X supp (j, U a (uj) := {uj k } k( - A 

A 

For a measurable set A C CI we denote by A A the cylinder set 

A* A := IT^IIaA) = {lj £ 0| 3c/ £ A such that U a (uj') = U a (lj)} 

The following observation plays a crucial role in the enhanced multiscale analysis. 

Lemma 4.8.1. For two disjoint cubes A, A' and two events A, B £ fL the induced 
events A* A and B* A , are independent. 

The following lemma allows one to turn a 'usual' Wegner estimate, as we have it 
discussed before, into a uniform Wegner estimate. It relies on the polynomial decay 
of the single site potential Q4.22)) , Let / be a compact interval, E £ I and e E]0, 1]. 
We denote by A(E,e,l) the event {u)\d(E,cr(Hl ) )) < e} and use the abbreviations 

n ; :=n Ai ,A* -.= A* Ai . 

Lemma 4.8.2. Under the above assumptions there exists a finite constant c, 
independent of lj G VL, l,r G N and e < 1 such that 

P{A(E, e, l)* l+r ] < P{A(E, e + cr^ m ~ d \ I)} 

Proof. By definition, for anwe ^l+ r there exists an u>' £ A such that 

Il2+ r a/ = ll; +r a; 

Thus, the coupling constants of ui and u)' with index k within the cube of size I + r 
coincide and we have for x £ A; 

\V u (x)-V u ,{x)\ < J2 \uk-u' k \u{x-k)<c' J2 \x ~ k\~™ < cr-( m -V 
\k\oo>l+r |feU>i+'' 

Therefore d(E, a(H l u ,)) < e implies d(E, cr(H l J) < e + cr-^ m - d \ which proves the 
lemma. □ 
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Let us have a look on the implications of the preceding lemma for a concrete 
example. Assume that the single site potential is bounded below on the unit cube 
around zero by k > 0. Then we have by Theorem l5.U.ll 

F{u\a(Hl) n [E - e, E + e] ^ 0} < C W (I) e l d 

for all E, e such that [E — e, E + e] is contained in the open interval /. This Wegner 
estimate implies its uniform analog 

(4.23) P({uj\a(Hl) n [E e, E] f 0}; +r ) < C W (I) (e + «■-<"*-">) l d 

for sufficiently large r > 0. In the application in the multiscalc analysis, both e 
and r are chosen as functions of I. From the estimate in (|4.23|) it might seem to be 
sufficient to choose m > d. This is also the minimal requirement to make the alloy 
type model with long range single site potentials well defined as a densely defined 
operator. However, for technical reasons, for the multiscale analysis to work one 
has to assume at least m > 2d. Under this assumption one can prove that the 
spectrum of H u is almost surely pure point near its bottom, and the corresponding 
eigenfunctions decay faster than any polynomial, 168, 329 . To obtain exponential 
decay of the eigenfunctions, one has to require m > 3d |169| . 

In the paper [329] by Zenk the above results have been extended to a model 
which incorporates random displacements of the single site potentials. 

5. Lipschitz continuity of the IDS 

In [1931 Kotani and Simon extended to continuum alloy type models certain 
arguments previously used for the derivation of Wegner's estimate for the discrete 
Anderson model. They treated only the case where the single site potential is the 
characteristic function of the unit cube, but Combes and Hislop showed in 54] that 
the same argument extends to non-negative single site potentials with uniform lower 
bound on the unit cube. There also some steps of the proof have been streamlined. 

One of the ideas in |193| is that in the same way as rank one perturbations are 
used for discrete Laplacians, positive perturbations may be used in the continuum 
case. This is related to the Aronszajn-Donoghue Theory [161 I17L 1181 I83j . See 
|45L 11921 12831 1282j for more background and references. This was essential, since 
a finite rank potential in the continuum may be a Dirac distribution, but not a 
function. 

Theorem 5.0.1. Let H u as in Definition \1.2.1\ and assume additionally that 
there exists an k > such that 

U > KX[_l/2,l/2] d 

Then for all E £ R there exists a constant Cw = Cw{E) such that for all I G N 
and all intervals I c] — oo, E] 

(5.1) E{Tr[P'(/)]}<CV \I\l d 

Remark 5.0.2. (a) It is sufficient to prove the theorem for the case k = 1. 
Since loqu = kujq k~ 1 u, the general case follows by rescaling the coupling constants 
and single site potentials. 

(b) The statement of the theorem remains true if one uses Neumann or periodic 
boundary conditions for H^, 
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(c) An explicit formula for the Wegner constant Cw is given in l|5.17jl . Since 
(|5~T|) is linear in the volume it follows \N(E 2 ) - N(Ei)\ < C w \E 2 - E x \. Thus, 
as we discussed already in § 13.11 the density of states n(E) := dN(E)/dE exists 
almost everywhere and is bounded by n(E) < Cw(E 2 ) for all E < E 2 . 

The next four paragraphs are devoted to the proof of Theorem I5.U.1I Up to 
some modifications we follow the line of argument in Section 4 of |54| . 

5.1. Partition of the trace into local contributions. In the present para- 
graph we derive preparatory estimates on 

(5.2) E{TrPi(/)} 

where we do not yet use the specific alloy-type structure of the potential. They 
have two aims. Firstly, to decompose the trace to contributions of unit cubes 
in A;. This will later facilitate the averaging procedure with respect to random 
parameters, whose effect on the potential is felt only locally. Secondly, it allows us 
to reduce the averaging of the trace of the spectral projection to the averaging of 
the quadratic form of the resolvent. The latter is technically easier to perform. 

Denote by A' and the Laplace operator on A; with Dirichlet, respectively 
Neumann boundary conditions. In § 11. 21 we saw that the potential V = V pcr + V u is 
infinitcsimally bounded with respect to — A and that the constants in the bound can 
be chosen uniformly in lu G fi. This implies that V is infinitcsimally form bounded 
with respect to any of the operators —A, —A' and — A l N with bounds uniform 
in w £ !1, 1 £ N and the choice of Dirichlet or Neumann boundary conditions. 
Consequently, there is a Co < oo such that for all oj S fl and I G N 

1(0,^)1 < ^,-A^)+C7 |H| 2 

which implies 

(5.3) (<t>,Hl<l>) > (0,™A^) - CII0II 2 > -Coll^l 2 

Thus + Co is a non-negative operator. 

Definition 5.1.1. A monotone decreasing, convex function r: [0, oo[-^]0, oo[ 
such that 

(5.4) C Tr := d,(r) := £ rfon))<oo 

will be called trace regularising. 

Throughout the rest of this section we denote by A the unit cube centred at 
zero. 

Remark 5.1.2. The bound i|5.4J) means that the operator r(-^A^) has finite 
trace. Namely, the eigenvalues of the Neumann Laplacian on the unit cube are 
given by 

d 

it 2 n 2 for all n £ Z d such that rij > 0, j = 1, . . . , d 

3=1 

cf. for instance |256j . page 266. By the spectral mapping theorem the eigenvalues 
of r(-±A£) are just rQ^E^nf). 
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Examples of functions r which are trace-regularising are the exponential func- 
tions r : x <— > e~ tx for t > 0. They have been used in |54l to implement the 
procedure outlined in this section. Another choice for r is a sufficiently high power 
of the resolvent x i— > (x + l)~ fc for k > d/2, which was used in |193| . That the 
operator x i— > (— ^Aj^ + l)~ fe is actually trace class can be inferred from |279| . 

The possibility to choose r from a large class of functions is of interest if one 
wants to give explicit upper bounds on the density of states. For instance, Section 
3.2 of |141] is devoted to deriving such explicit upper estimates. However, there, 
following |54 |. the function r(x) — e~ tx is used. Due to this choice, the upper 
bound on the density of states is exponentially growing in the energy. This can 
be improved to a merely polynomial growing bound. Furthermore, if one studies 
coupling constants which may take on arbitrarily negative values, the choice of r 
will determine which moment conditions one has to impose on the negative part of 
ujo, see also § 15.61 

Proposition 5.1.3. With C as in (15.31) 

E {TrPi(I)} < r(E 2 + Co)- 1 C Tr (r) ]T || E {x^0OXi}|| 

jeAi 

Proof. Since =■ is well-defined and bounded on the compact interval / := 
[Ei, E2], we have 

Tr [P< (/)] = Tr [r{H l u + Co)" 1 P< (/) r(H l u + C„)] 

Furthermore, by spectral calculus and since for positive operators A, B we have 
Tr(AB) < ||A|| Tr(B), the above line is bounded by 

r(E 2 + Ca)- 1 TT[Pl(I)r(Hl + C )} 
According to the direct sum decomposition 

i 2 (A,)=0i 2 (A + j) 

we consider the Laplace operators —A J,N on L 2 (A + j) with Neumann boundary 
conditions. Dirichlet-Neumann bracketing implies 

(5-5) Hi + C > - l -A l N > -\ A^ N =: ®H 

For a normalised eigenfunction 4> of corresponding to the eigenvalue A we have 
by the spectral mapping theorem 

(5.6) (cf>,r(Hl + Co)<f) =r(A + C ) - r((4>, (H l u + C )<f}) < r((<p,®H<t>)) 

Applying Jensen's inequality to the spectral measure of @H we estimate (|5.6|) from 
above by (</>, r{®H)4>). Let 4> n > n G N be an orthonormal basis of eigenvectors of 
with corresponding eigenvalues A„, n £ N. We apply the above estimates to the 
trace 

Tr [Pl{I)r{Hl + C Q )j < £ (cf> n , r(H l u + C )<M 

< (<Pn,r((5H)<p n ) <Tv[Pi(I)r((BH)] 
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For the next step we write down the trace with respect to different basis. For each 
j £ A let {^£1 n £ N} be an orthonormal basis of L 2 (A+j), then {ip J n \ n £ N, j £ A} 
is an orthonormal basis of L 2 (Ai). Since r(®H)ip^ — Xj r {~\^' N )Xj' i l4i ^ follows 
for the trace 

= E E^>'Xi-P 1 i(/)XiK-|A i ' JV )Xi^>) 

jgA, n£N 

jeAi 

Thus we have decomposed the trace to contributions from each unit cube. We 
summarize the estimates so far: 

TrP< (J) < r(E 2 + Co)- 1 ]T Tr [ Xj P^(I) Xj r(-±A^) Xj ] 

Since I is a bounded interval and V peI + Kj is an infinitesimally small perturbation 
of —A' independently of ui, it follows that the dimension of P^(I)L 2 (Ai) is bounded 
by a constant C3 uniformly in ui. Thus 

Tr [XjP l M)xA-W' N )Xi\ < C 3 r(0) for all w e fi 

is an upper bound by an integrable majorant and wc are able to interchange the 
trace and the expectation by Lebesgue's theorem on dominated convergence 

E {Tr [^P*(I) Xj -r(-fA^) Xj -]} = Tr [E { X] Pl{I) X3 r{-h^ N )Xj}] 

= Tr[E { Xj Pl{I) Xj } X A-\^' N )Xo] < ||E {XjPLWXi}\\ Tr [r(-|A»- w )] 

By assumption, r is trace regularising, so the trace in the last line is finite. □ 

5.2. Spectral averaging of resolvents. Now we consider how resolvents are 
averaged when integrated over a random parameter. Together with the partition 
result in the previous paragraph § 15.11 this will enable us to complete in § 15.41 the 
proof of Theorem l5.0.1l 

Apart from this, the spectral averaging result bears in itself a meaning. Con- 
sider a nonnegative operator H with discrete spectrum. Its resolvent R{E) = (H — 
E)- 1 has singularities at the eigenvalues of H which are of the form (A n — -E) -1 , A n £ 
<j{H) and thus are not integrable over the energy axis. In other words, for a general 
vector (j> the function E 1— > (</>, R(E)tfi) will not have a convergent integral. Now, 
if H = H\ depends on a random parameter A, we might hope that the averaged 
resolvent E 1— > J dP(X)(<f>,Rx(E)<f)) will be integrable. This would mean that the 
singularities of the resolvent have been smeared out sufficiently by the integral over 
A. The lemma in this paragraph shows that this is actually the case for operators 
which depend in a specific way on the random parameter. 

Consider the following operators on a Hilbert space 7i. Let H be a selfadjoint 
operator, W symmetric and infinitesimally bounded with respect to H, and J non- 
negative with J 2 < W. Choose two parameters 

z £ C_ := {z £ C| Imz < 0} 

(£C^:={(£C\ ImC>0} 
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and set 

(5.7) H(Q :=H + (W, K((, z) := J(H(Q - z)~ 1 J 

The following lemma is a slight generalisation of Lemma 4.1 in |54| . 

Lemma 5.2.1. For all z G C_, all t > and any normalised </> G TL we have 

d( 



(5.8) 



l + *c 2 



< 7T 



Proof. By Pythagoras we have \{4>, (A + iB)4>)\ 2 = |(0,v40)| 2 + |(0,Bc/>)| 2 
for any two selfadjoint operators A, B. Thus the norm of K (£, z) is bounded by 
| Imz| _1 1| J|| 2 . On the other hand, the equation 



implies 
(5.9) 



■ hnK((, z) = B[(H(C) - zT 1 [(Im QW - Im z] (H(C) 



i^c^Hiimcr 1 



B 



Here we used that W(H(() — z)^ 1 is a bounded operator. Now observe that for all 
z G C_ the function £ i— > K{Q, z) is holomorphic and bounded on C+. The residue 
theorem, integration over a closed curve in C and the bounds on K imply 



(5.10) 



(0,if(C,z)0) 



Vt 



\\K(i/Vt,z 



Together with l|5.9|) . this completes the proof. 



□ 



Remark 5.2.2. The lemma shows that for the particular family of operators 
-ff(C) in H5.7fl . where £ is a random variable with measure pp^a ; the C-averaged 
resolvents are indeed integrable with respect to the energy. Thus the singularities 
of the resolvent have been smeared out. 

5.3. Stone's formula and spectral averaging of projections. Stone's 
formula allows one to express the spectral projection in terms of the resolvent. 
This is handy because the resolvent has some nice analytic properties. In our case 
we use Stone's formula to derive the analog of l|5.8|) for spectral projections. 

A sequence of bounded operators A n , n € N on the Hilbert space H converges 
strongly (or in strong topology) to A if for every tf> £ Ji 



lim \\A(j> - A T: 







Lemma 5.3.1 (Stone's formula). Let H be a selfadjoint operator with spectral 
family denoted by P(-). Then the following limit holds in the strong topology 



lim / \{H — E 

s\o2TTiJ El L 



iS)- 1 - (H-E + iS)- 1 ] dE 

1 r 



P([E 1 ,E 2 \)+P(]E U E 2 [) 
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PROOF. The function 

If x — Ei x — Ei 
(5.11) js{x) := — arctan arctan — 



7T 



f *\(x-E- id)- 1 -(x-E + iS)- 1 } dE 
2tt« J El 

1 f® 2 

Im / (x- E + iSy 1 dE 

n Je 1 



converges for 5 \ to 

-^(X[Ei,E 2 ] + X]E 1 ,E 2 [)- 

Now one applies the spectral theorem to fs(H). □ 

More details on Stone's formula can be found in |258| . or |326j where the 
spectral calculus is actually introduced in this way in Section 7.3. 

Now let H (£) be as in the last paragraph and P(£, /) the corresponding spectral 
projection onto an interval I. For a normalised vector tp in TL denote V (£) :— 

JP(C, I)Jip) ■ The next lemma contains a spectral averaging estimate for V . 

Lemma 5.3.2. Let p € Z°°(R) n L 1 (M). Then 
(5-12) /" p(C)^(C)dC< IHIoo|/| 



While Combes and Hislop 54 considered only compactly supported p, it was 
first observed in jlU7j that densities with non-compact support can be treated. 
There this extension was necessary to derive estimates for Gaussian random poten- 
tials. 

Proof. We first consider the special density j^^i and an open interval I. By 
Stone's formula 

(5 ' 13) L IT? v (c) * - L TW & \ Im / dE ^ K{c > E ~ ^ 

Note that j^§^2 is a finite Borel measure on M and that (|5.11|l implies that |/$(-)|, 
and hence \\fs(H (C))|| 5 is bounded by one. Thus we may apply the dominated 
convergence theorem to interchange the limit and the integration, and bound (|5.13|) 

by 

(5.14) - lim 

The last inequality follows from Lemma [5.2.11 This implies for all p £ L°° with 
compact support: 

P (0r(C)d( < sup [ P (0(i + iC 2 )] I ^M*<K 



< \I\ 



suppp 

< sup [p(C)(l + ^C 2 )] |/ 

suppp 

-> \\pWoo\I\ fort^O 



i + K 2 



Finally, assume only that p <E L°° n L . Set p v = px{x\ \x\<y} and decompose 
p = p y + p y . For y — > oo, p y tends to zero pointwise. Since V is bounded by one, 
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p € L (R, d£) is a y- uniform majorant for p y V and we may apply the dominated 
convergence theorem to conclude 

(5.15) f p(QP(Od(= lim f py(OT(()dt< ||H|oo|/| 

Jm y^°° Jr 

If / is not open, we write it as an intersection of open, decreasing sets and use 
monotone convergence to conclude (I5.12|) . □ 

5.4. Completion of the proof of Theorem 15.0.11 The results on the lo- 
calisation of the trace to unit cubes and spectral averaging of projections allow us 
to assemble the proof of Theorem 15. 0.11 

To estimate the operator norm appearing in Proposition 15 . 1 .31 we may as well 
bound the corresponding quadratic form since 

(5.16) || E {xjPL(I)Xj}\\ = sup (0, E { Xj P l u (I) Xj } 4>) 

11011=1 

Now one can apply Fubini's Theorem and Lemma 15.3.21 with the choice p = f, 
H = i 2 (A z ), J = X j, H = H Q + Efc 6 A\j u k u(- - k), ( = tOj and W = u(x- j): 

(0,E { Xj Pi (i)xj} 4>) < \\f\\oo\i\ 

This bound is j-indcpendent and thus yields 

(5.17) E {TrP< (/)} < r(E 2 + Co)' 1 C^(r) ||/|U \I\ \A t \ 

□ 

Now it becomes clear why we introduced the operator r(ii^ + Co) _1 r(i?^ + 
Co) = Id in Proposition 15. 1 .31 without this regularisation of the trace we could 
have estimated 

E{Tr Xj PUl)Xj} < const\I\\Ai\ 
However, this would lead to a Wegner estimate with quadratic volume bound. 

The role played by r resembles the one of the function g in paragraph 14.21 and 
Appendix 1X1 

5.5. Single site potentials with changing sign. In § 14.71 we saw an ex- 
tension of the Wegner-Kirsch approach to single site potentials of changing sign. 
The Kotani-Simon-Combes-Hislop proof of Wegner's estimate also allows such a 
generalisation [3151 1317) . which we present in this section. Its main shortcoming 
in comparison to the results in § 14.71 is that it is restricted to single site potentials 
which have a (generalised) step function form. On the other hand, it is valid not 
only at spectral boundaries, but on the whole energy axis. Furthermore, it yields 
the existence of the density of states as a function and upper bounds on it. 

Theorem 5.5.1. Let L p c (M. d ) 3 w > nX[o,i] d k > and p(d) be as in 

(|1.1() . Let r C Z d be finite, the convolution vector a = (a^ker £ satisfy 
a* := Efc^o l a fe| < l^o I j an d the single site potential be of generalised step function 
form: 

(5.18) u(x) — otk w(x — k). 

feer 

Assume that the density satisfies f € W^' (M)- Then for all E £ R there exists a 
constant C'w = CV(-E) such that 

(5.19) E {TtPI{L)} <C W \I\l d , for all I and 7 c] — oo, E] 
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The theorem implies that the DOS, the derivative of the IDS, exists for a.e. E 
and is locally uniformly bounded: dN(E)/dE < C{E\) for all E < Ei. 

Proof. For simplicity we assume w — X[o,i] d - To estimate E {((f), XjPLi^Xj^)} 
for any normalised 4> £ L 2 (Ai) we introduce a transformation of coordinates on the 
probability space f2. 

For each cube A = A ; denote A+ := {A - 7] A S A, 7 e T} and L = #A + . The 
operator depends only on the truncated random vector (oJk)keA+ G On 
such vectors acts a multi-level Toeplitz matrix A\ := {a>j-k}j,keA+ induced by the 
convolution vector a. The transformation has an inverse -Ba = {^fc,j}fe,ieA+ = 
which is bounded in the row-sum norm ||Ba|| < i- a * • Note that the bound is 
uniform in A;. We drop now the subscript A and denote with 77 := Auj the vector 
of the transformed random coordinates. They have the common density 

(5.20) k(r}) = I dct B\ F(A^ 1 ?]) 

where F(u) = IlfeeA+ i s the original density of the w&. We calculate the 

potential written as a function of r\ (and x e A): 

V u {x) = V Bv (x) =^2vjXj(x) 

In the new representation of the potential the single site potentials are non-negative, 
so we can make use of the spectral averaging formula in Lemma 15.3.21 

(5.21) f d Vj k( V ) s( V ) < \I\ Bup|%)|, where 5(77) := (</>,XjPL,(I)Xj4>) 

JR Vj 

Fubini, (|5.21l) . and the fundamental theorem of calculus give 

(5.22) f drjk(rj)s(r])<\I\ I drj ±j sup |fc(r?)| < |/| / drj \{djk)(ri)\ 

Here is an abbreviation for \r\h\ k e A + \ j}. The last integral equals 

I det A| J du\{djk)(Auj)\ 

which is bounded by H/'Hl 1 J2keA+ \bk,j\- The proof of the theorem is finished by 
the estimate 

(5-23) E{(^ Xj Pi(I) X ^)}<\I\ \\f\\ L 4B\\ 

and Proposition 15.1.31 □ 

One drawback of Theorem l5.5.1l is the requirement of the weak differentiability 
of /. This excludes in particular the uniform distribution on an interval. However, 
in a joint work |191| with Kostrykin we have proven: 

PROPOSITION 5.5.2. Let the assumptions of Theorem \5.5.1\ be satisfied with the 
only difference that f is the uniform density on an interval. Let T C {k £ Z d | fcj > 
V i = 1,. . . Then l|5.19|l holds true. 

In fact, it turns out that the proof of Theorem 15. 5. II can be extended to density 
functions of finite total variation. This covers in particular linear combinations of 
functions in W^' 1 and (finite) step functions. More precisely: 
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Proposition 5.5.3. Let the assumptions of Theorem \5.5.1\ be satisfied, but 
require f merely to have finite total variation ||/||var < oo. Then the Wegner 
estimate l|5.19|l holds. 

The difference to Theorem 15.5. II is that the constant Cw now depends on 
ll/llvar instead of cf. (|5.23|1 . The result in the Proposition 15.5.31 is proven 

in |190j . Moreover, there we discuss how the condition J2k^o \ ak \ < l a o| can be 
relaxed using the theory of Toeplitz matrices. Let us give an example in the one 
dimensional case. 

The symbol of the Toeplitz matrix A is the function 

s A : T^C, s A (e l9 ) = £>,e iie , 8 e]-ir,n] 

jez 

Since we assume that only finitely many components of a are different from zero, sa 
is actually a trigonometric polynomial, and thus uniformly continuous and bounded. 
Invertibility criteria for A\ as well as bounds on B\ = A^ 1 and B = A^ 1 may be 
established by studying the symbol sa- Consider the case that the symbol sa has 
no zeros and the winding number of sa with respect to G C vanishes. A theorem 
of Baxter EU, see also (T551 Thm. III. 2.1], states that 

(5.24) sup(||B||,||Ba,||) < const. < oo 

In this case we have again: 

Proposition 5.5.4. Let d = 1 and the assumptions of Theorem 1 5. 5. 1\ be satis- 
fied, but require for the convolution vector a — (a^ker merely that the symbol sa 
of the associated Toeplitz matrix has no zeros. Then the Wegner estimate (|5.19|) 
holds. 

Remark 5.5.5 (Anderson model). For the discrete Anderson model h u = ho + 
Vuj there is a result analogous to Theorem 15.5.11 Here ho is the finite difference 
Laplacian on l 2 (Z d ) and (V^VOC 71 ) = V^(n)'0(:re), Vn E Z rf , a multiplication operator 
as in the continuum case. This is not surprising, since the arguments in § 15.21 and 
t|5.3l rely only on abstract functional analysis. If fact, as we mentioned earlier, 
Kotani and Simon were motivated in their treatment |193| of the alloy type model 
by its discrete counterpart. Moreover, since on l 2 (Z d ) the trace can be expressed 
using the canonical basis as 

the use of a trace regularising function is not necessary. Here denotes the 
spectral projection of the truncation h 1 ^ of the Anderson model h u . More precisely, 
hi is the finite matrix {{Sj, h^d}.)} j ke ^. 

Note that in the discrete case Xj is j us t Sj. Under the assumptions of Theorem 
15.5.11 on the coupling constants {ujj }j and the single site potential u we have the 
following Wegner estimate for the Anderson model: 

(5.25) E {Tr[Pi(/)} < JLOg |/| |A,| 
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Remark 5.5.6. Theorem 15 . 5 . II can also be understood as a Wegner estimate 
for the alloy type potential 



where the coupling constants {f]j}j are not any more independent, but correlated 
satisfying certain conditions. See § 4.2 in |318| for a precise formulation. Wegner 
estimates for correlated coupling constants can also be found in |59| (cf. |141| . too). 

The use of the common density F, respectively fc, in the proof of Theorem 
15.5.11 is conceptually new. One could try to use conditional densities instead by 
considering the indefinite potential V u in its representation Vbt? as an alloy type 
potential with dependent coupling constants. However, this would require to have 
uniform upper bounds on the conditional densities, cf. |59L 1141] . They do not 
seem to be easy to establish for the model considered in this paragraph, and in fact 
sometimes fail to hold as can be seen in the following example. 

Example 5.5.7. It is sufficient to consider only one space dimension d = 1. Let 
the density function be / = X[o,i] an d the single site potential u = X[o,i] ~ a X[i,2] 
with —a e] — 1, 0[. To this model the results of Propositions 15 . 5 .21 and 15331 apply. 

The restriction of to the interval ] — 1/2, 1 — l/2[ of length I depends only 
on the coupling constants ujj with indices j £ {— 1, . . . , Z — 1}=: A + . They are 
transformed by the Toeplitz matrix A into new random variables rjj, j G {— 1, . . . , I— 
1}, as in the proof of Theorem 15.5.11 Here the convolution vector is given by 
ao = 1, ct\ = —a. 

The conditional density Pj(i]) — pArj) of the variable rjj with respect to the 
remaining coupling constants = (rjk)keA+\j m A + is given by pj(rj) — ■ 
Here gjirj) — J k(r))dr]j denotes the marginal density. The question is whether 
supj p(rj) is finite. 

One calculates the common density to be k{rf) = nl=-i X[o,i](X^=-i ot k ~ u rj v ). 
For r)j + i e [0, 1], ?7fe = 0, Vfc 7^ j + 1 we have 




feez d 



k(rj) 



k=j+i 



The marginal density 





has for rjj + i G [0, 1],% = 0,Vfc ^ + 1} the upper bound 




Particularly, gjijq) \ for Jjj+i Z' 1 and thus 

sup pj (rf) = oo 
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Therefore, proofs of a Wegner estimate which require the conditional density to be 
bounded cannot be applied to this alloy type potential. See also § 4.3 of 318 for 
another example. 

5.6. Unbounded coupling constants and magnetic fields. Motivated by 
certain physical applications, e.g. the study of the quantum hall effect (see for 
instance |26l 12341 11111 12711 11531 I92| ). it is desirable to extend the results on 
the continuity of the IDS to include Hamiltonians with magnetic fields. This is, for 
instance, done in the papers |SQJ I134L I141j . 

We discuss here the results on alloy type potentials obtained in |141| by Hupfer, 
Leschke, Miiller, and Warzel, since they are build on the method presented in the 
preceding §§ I5.1H5.4I Moreover, their result allows the coupling constants to be 
unbounded, as long as very negative fluctuations are exponentially rare. Actually, 
the primary interest of their research are Hamiltonians with Gaussian random po- 
tentials, so they need to cope with unbounded fluctuations of the potential. The 
proof is based on earlier techniques from |1(J7) — which in turn use |54| — and 
Dirichlet-Neumann bracketing for magnetic Schrodinger operators, as discussed in 
Appendix A of |141| . The results concerning alloy type potentials are summarised 
in § 4.1 of their paper, which we review shortly. 

Let A: R d -> R d be a measurable vector potential with the property \A\ G 
L 1 1 oc (R d ). Denote with H Q the selfadjoint closure of Yji=j(i®j + Aj) 2 defined on 
smooth functions with compact support. The alloy type Schrodinger operator = 
Ho + V u now incorporates a magnetic field. In |141| it is proven that Thcorem l5.U.ll 
essentially remains true if the magnetic field is included. (Their conditions on the 
single site potential are slightly different.) 

Moreover, two cases are discussed, where the coupling constants are unbounded 
random variables, and the theorem still remains true. In the first one, it is as- 
sumed that ujq is non-negative and a certain moment condition is satisfied, roughly 
E{c^+ 2 } < oo. 

The second one concerns the case where luq is distributed according to the 
Laplace distribution: P{ujq G /} = - Jj dx e^ x ^ a . Note that the probability that ujq 
assumes very negative values is exponentially small. In fact, in |141| it is noted, 
that this is a necessary requirement for their techniques to work. The reason for 
this is that they use r(x) — e~ tx as the trace regularising function, cf. Definition 
15.1.11 A different choice of r would allow for more general distributions unbounded 
from below. 

We conclude this section by listing further literature on random Hamiltonians 
with magnetic fields. Works treating the regularity of the IDS of random Schro- 
dinger operators with magnetic field include |321L 1221 13231 11451 1134j , while the 
question of the (in)dependence of the IDS on boundary conditions for these models 
has been treated in |232l 1821 11411 1139] . A related problem is the analysis of the 
semigroup kernels of magnetic operators |38l I39 |. In 166, 1253| the behaviour of 
the IDS in a strong magnetic field is identified. 

The asymptotic behaviour of the IDS near the boundaries of the spectrum in 
the presence of random magnetic fields was the object of study of the articles |230L 
12311 1233) which prove high energy and Lifshitz asymptotics for certain models. 
The high energy asymptotics has been analysed already in |221L 1312) . 
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For Schrodinger operators with constant magnetic field and random potential 
generated by a Poissonian process the different possible behaviours of the IDS at 
the bottom of the spectrum are analysed in [3711551 H151 ITIil 1351 15211 IT35] . 

The analysis of Landau Hamiltonians in the single band approximation is done 
in [841 1851 12461 1247L I145| . Examples of localisation proofs which allow for mag- 
netic fields can be found in [84l f55ll85ll86l l322l l22l ITT81 RTfl H09l 12471 16]. 

Appendix A. Properties of the spectral shift function 

For a exposition of the theory of the spectral shift function (SSF) see |33l 1282] 
or the last chapter of |328j . We review here the relevant facts in our context. 

For two trace class operators A, B the SSF £(-,A,B) may be defined by the 
formula 

(A.l) Tr(f(A) - f(B)) = J /'(A) £(A, A, B) dX 

for functions / £ . Actually, it holds for more general functions, too. By 
Theorem 8.3.3 in 328 , it is sufficient to assume / S C 1 (IR) and that /' is the 
Fourier transform of a finite complex measure. One can define the SSF also via the 
perturbation determinant from scattering theory 

£(A, A, B) := - lim argdet[l + {A - B)(B - A - ie)" 1 )} 

In this case the equality IjA.ljl is called Krein trace formula. The SSF can be 
bounded in terms of the properties of A — B, namely 

(A.2) ||£(A,A,B)|| 1 <|| J 4-.B|U 

Here J\ denotes the ideal of trace class operators and || • \\j 1 the trace norm. On 
the other hand, if A — B is finite rank 

(A.3) ||£(A,A,£)|| 0O <mnk(A-B) 

Since we have an estimate on £ in the L 1 and L°°-norms, it is natural to ask whether 
an estimate for the L p -norm, p g]1, oo[, may be derived. This indeed turns out to 
be true and can be understood as an interpolation result, cf. the proof of Theorem 
2.1 in 160) . 

To formulate this bound we have to introduce ideals of 'better than trace class' 
operators. For a compact operator C denote by fi n {C),n £ N its singular values, 
in non-increasing order. If C is trace class, the sum of the singular values is finite 
and equals ||C|| j x . We denote by Jp the class of compact operators such that 

(A.4) \\C\\j, :=( £ Mn(Cf) <oc 

neti 

The theory of such operators is classical for f3 > 1. However, since we want to 
interpolate between l|A.2(l and (|A. 3|) , we need to consider operators whose singular 
values converge not slower, but faster than a ^-sequence to zero. This leads us to 
consider operators such that ||C||jg is finite, for (3 smaller than one. In particular, 
all such operators are trace class, which explains why they are sometimes called 
super-trace class. It follows that the SSF may be defined for such operators. They 
have been studied in |126ll5Ull31) . while their relevance in the present context was 
recognised in |6U| . 
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Form these sources we infer the following properties of Jp . Since for any com- 
pact operator A and bounded B the singular values of the products obey 

(A.5) n n (AB) < \\B\\ n n (A) and (i n (BA) < \\B\\ fi n (A) 

the set Jp is an two-sided ideal in the algebra of bounded operators for all (3 > 0. 
For (3 > 1 the functional A i— > ||A|| is a norm, which is not true for (3 < 1. More 
precisely, in this case we have only 

\\A + B\\%<\\A\\% + \\B\\% 

This property implies that || • is a quasi-norm and that 

distp{A,B) = \\A-Bf Jf} 

is a well defined metric on Jp. The pair (J^, dist^) forms a complete, separable 
linear metric space, in which the finite rank operators form a dense subset. 
In 6CQ the following L p -bound on the SSF was proven. 

Theorem A.l. Letp> 1 andA,B be self adjoint operators whose difference is 
in Jp where (3 — 1/p. Then the spectral shift function £(-,A,B) is in L P (R) and 

(A.6) U(;A,B)\\ L ,<\\A-Bf Je 

This estimate is sufficient for our purposes. There exists a sharp version proven 
by Hundertmark and Simon in 140 . It is used in the result described in Remark 

Theorem A. 2. LetF: [0,oo[— > [0, oo[ be a convex function such thatF(0) = 0. 
Let A, B be bounded and C a non-negative compact operator such that for all N G N 



(A.7) $> n (|A-B|)< $>n(C) 

n—N n—N 

Then 

F{\i{\A,B)\)dX < f F(\t[(\,C,0)\)dX = i F ( n ) - F ( n - !)] ^{O 

Condition (|A.7|I is in particular satisfied if \A — B\ < C. Of course, to apply 
Theorem lA.il we need a criterion for the operators which arise in our situation to 
be in Jp for (3 — 1/p < 1. So, let's have a closer look at the application of the 
theory of the SSF to Schrodinger operators. 

Since we are studying Schrodinger operators, we cannot expect to deal with 
trace class perturbations. However, the theory extends to operator pairs such that 
the difference of a sufficiently high power of their resolvents is trace class. More 
precisely, assume that H + u, H is a pair of lower bounded operators such that 

(A.8) H + u > C , H > C and g(H + u) - g(H) e J l 

where g(x) — (x — Co + l)~ k for some sufficiently large k > 0. Following 60], we 
denote g(H + u) — g(H) by V e s. This is the 'effective' perturbation, although it is 
obviously not a multiplication operator. One defines the SSF of the pair H + u, H 



(A.9) Z(\,H + u,H):=-Z(g(\),g(H + u),g(H)) for A > 
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and £ = otherwise. This definition of £ is independent of the choice of the 
exponent k > in g. By Theorem 8.9.1 in |328| the trace formula (|A.1J| holds if 
/ e C 2 (M.) and /' has compact support. This conditions are clearly satisfied by the 
switch function p we use in § 14.11 and § 14.21 

The purpose of the theorem we are heading to now is twofold: firstly, to estab- 
lish that V e ff is in Jp for suitable (3 < 1. Thus, we will be able to apply Theorem 
IA.1I Secondly, to control the upper bound ||V^fF|| j appearing in (|A.B|) . 

It is well known that operators which may be formally written as f(x)g(—iW) 
are in the Hilbert-Schmidt class if /, g E L 2 (M. d ). The product of such two operators 
is trace class. Extending this idea, we want to show for certain operators that they 
are in some (super-trace) ideal Jp , (3 < 1 by writing them as a product of sufficiently 
many operators of the type /(x)</(— zV). For this purpose it is useful to note that 
the Holder inequality extends also to the case of exponents smaller than one: let 
a L : N — > C, i = 1, . . . , N be such that |o,i(n)| Pi is summable, where pi > for all 
i = 1, . . . , N, and set £ := J2i=i Then the pointwise product fli=i a » ^ s m ^ r (N) 
and 

N i N 

inHi ^nii a «!ip.. 

i=l T i=l 

By applying this to the sequence of singular values of compact operators, we obtain 
the following 

Lemma A. 3. Let A{ G J Pi {or i — 1,..,,N, then Jli=i * s * n Jr where 

i 1 J 

- := ) — and 

N N 

(a.io) llrHL ^nii^ik 

i=l r i=l 

See also |31| Corollary 11.11. The following result is taken from 232 , cf. Lemma 
11 and its proof. For IgNwe abbreviate A = A;. 

Lemma A.4. Let q > d/2,q e 2N, / € Lf(R d ) and C € R be such that 
Vpcr + Vlo > Co /or aZZ w. TTien iTie operator product / (iT*, — Co + l) -1 is zn ifte 
ideal J q and 

H/XA (ffi - Co + 1)-%, < \\f(-A + l) _1 ||,/ g < C{q) 11/11, 

PROOF. There exist a bounded extension operator £ : W 2 - 2 (ki) — > VF 2 ' 2 (R d ) 
and its norm is independent of Z e N, cf. Section IV. 3. 2 in 293]. Thus we have 

fXA{Hl - Co + I)" 1 = f X A(H u - Co + l)- 1 ^ - Co + l)£(Hl - C* + l)" 1 

Since £ is an extension operator, xa£ is the identity on W 2 ' 2 (Ai). By the ideal 
property of J q and the boundedness of £ we have 

\\fXK (Hl-Co + l)- 1 ^ 

< \\f(H u - C + l)-% q \\(H U - Co + l)£(Hl - Co + l)" 1 ]! 

< const. \\f(H u - C + !)-%„ 
By the Kato-Simon |151L I277| or diamagnetic inequality we have 
(A.ll) \f(H u - Co + < /(-A + 1)-Vl 
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for ip E L 2 (M. d ), ci. proof of Theorem 3.3 in |139j . For the 'free' case we know by 
Theorem 4.1 in j2?8] / (-A + l)" 1 € J q and 

\\f(- A + \\f\\ g \\h\\ g 

where h(x) = (x 2 + 1) . The pointwise inequality (|A.11(I implies by Theorem 2.13 
in [2H] / (H u - Co + l)" 1 € J q and 

||/(^-C + l)- 1 ||.7 g <||/(-A + l)- 1 ||. /5 

□ 

Remark A. 5. The result remains true if we consider Neumann boundary con- 
ditions instead of Dirichlet ones, and if we include a bounded magnetic vector 
potential in the background operator Hq, see Lemma 10 in 232 . This fact re- 
lies on the existence of an appropriate extension operator which takes (magnetic) 
Sobolev functions on A to (magnetic) Sobolev functions on M d . 

The following result establishes that V c g = g(Hl+u + u) — g(iJ^+W-w) is indeed 
super-trace class and that its quasi-norm HV^ffHjp can be bounded independently 
of the cube A. 

Theorem A. 6. Let Hi = -A + V and H 2 = Hi + u, with ±C < V,u S 
Lf ocunii (R d ), where p is as in (11.1(1 . Denote by H[,H 2 the corresponding Dirichlet 
restrictions to the cube A;. Assume 

|| V|| p ,unif,loc < Ci, and ||lt||p,unif,loc < C 2 

and that the support of u is contained in the ball Br(x). For any (3 > choose 
k e N with k > 4±^. For g(x) = (x - C + l)- k set V eS := g(H 2 ) - g{H{). 

Then V e & € Jp and ||T4ff|| j. is bounded by a constant which is independent of 
A; and x and depends on V and u only trough Co, Ci, C 2 and R. 

PROOF. Choose a function / e C^°(B 2 r(x)) such that / = 1 on Br(x). An 
iteration of the resolvent formula yields 

k-l 

V cS = -J2 ( ff 2 - Co + l)-( fc - m > u (Hi - Co + l)-( m+1 > 
m=0 
k-l 

(A.12) = - [f k ~ m {H 2 - C + l)-( fc -™)] * u [f m+1 {Hi - Co + l)-( m+1 )] 

m=0 

By Lemma 9 and Appendix A in 232 we have the following representation 

N v 

r (Hi - Co + \)- v = n fa ^ - c " + i y lB *i 
i=i j=i 

and analogously for H 2 . Here TV = N(y) is an integer which depends only on v 1 
the functions /y € C^°(B 2 r(x)) are linear combinations of derivatives of /, and 
Bij are bounded operators with norms independent of I. Set 

C q ,f, B = C(q) max{||B«|| ||/y L . | i = 1, . . . , N, j = 1, . . . , v} 



INTEGRATED DENSITY OF STATES AND WEGNER ESTIMATES 



167 



where g£]f,f + 2] and C(q) are as in Lemma fA. 41 By the same lemma, the ideal 
property of Jp and the Holder inequality (| A. 1 0|> we have 

(A.13) 

N v 

wr (H t -c + i)-"n* Jt < E II n/« - c ° + ^"'IIS, ii^ii* ^ NC h,B 

i=l 3=1 

ift = q/v< 1. For t > 1, || • ||j t is even a norm and a similar estimate holds. Since 
it is infinitesimally bounded with respect to the Laplacian 

C u := max{||«/« (#i - Cq + l) -1 ^!!} 
is finite. Thus, in analogy to (|A.13|I 

\\ur +i (ifi - Co + i)- y+i ns f <»c: c* /iB 

From the preceding we see that l|A.12|) factorises V D s as a product of bounded 
operators and k operators which are in J q . All the involved operator and super-trace 
class norms can be bounded independently of A; and x. Using the ideal property 
and Holder's inequality we see that V c g is in J r for all r > q/k and ||14ff||j r is 
bounded by a constant which is independent of A; and x. The way we choose k 
and q makes it possible to take r — (3. □ 

Remark A. 7 (Properties and relevance of the SSF). An exposition of the role 
played by the SSF in scattering theory can be found in |328j . The SSF has proven 
useful in the study of random operators, particularly in problems related to surface 
models, e.g. the definition of the density of surface states |48l 1491 11881 1185) . 

Various of its properties are discussed in the literature: monotonicity and con- 
cavity |117L 11211 1187| . the asymptotic behaviour in the large coupling constant 
[2im l2l)5l [232] and semiclassical limit and some other bounds [2391 [2301 

1215]. 
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— regularising, 153, 155, 162 

unique continuation property, 131, 133, 149 

van Hove property, 115, 124 
volume 

— comparison theorem, 121 

— density, 109 

von Neumann algebra, 112 

Wegner constant, 130 

Wegner estimate, 129, 134, 135 

Weyl asymptotics, 127, 134, 136, 141 
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